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1. 


Introduction 


This  volume  of  the  report  serves  as  a  user’s  guide  to  the  Illinois  Institute  of  Technology 
(ITT)  Signal  Processing  Worksystem™  (SPW)  simulations.  These  simulations  consist  of  nine 
modular  SPW  systems  and  fifteen  Memoryless  Nonlinear  Processor  (MNP)  blocks  for  use  with 
these  systems.  A  solid  understanding  of  the  concepts  discussed  in  Volume  /  of  this  report  is 
essential  for  using  these  systems  and  interpreting  the  results. 


NOTICE  TO  THE  USER: 


The  most  important  system  parameter  relationships  in  this  report  are 
summarized  in  three  text  boxes  like  this  one.  It  is  imperative  that  these 
fundamental  constraints  are  maintained  at  ail  times  during  simulation 
of  the  systems.  There  are  no  provisions  in  the  systems  to  automatically 
verify  these  constraints;  it  is  entirely  the  responsibility  of  the  user.  If 
these  constraints  are  violated,  invalid  results  may  be  produced. 


2.  The  HT  SPW  MNPs 

In  the  current  and  past  research  efforts  five  types  of  MNPs  have  been  implemented  in 
the  SPW  platform:  Histogram  [Illi93a],  Equiprobable  Bin  Histogram  (EBH)  £Illi93b] ,  Fourier 
Series  Approximation  (FSA)  [Illi93b],  Continuous  Polynomial  Approximation  (CPA)  [Grim93], 
and  M  Interval  Polynomial  Approximation  (MIPA)  [UH91].  Table  1  summarizes  the  filenames 
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for  the  MNPs'  and  the  linear  receiver.  These  filenames  are  for  reference  only  since  they  may 
be  copied  directly  from  the  SPW  User  Palette  [Comd91]  into  the  desired  system. 


Filenames  of  MNP  Blocks 


NONLINEAR  RECEIVERS  WITH  CORRELATOR: 


Librarv/Filename 

HISTO/SERIALNONLIN 

HISTO/S_ENONLIN 

FOURIER/SERIAL_NONLIN 

CPA/S_CDF2_NONUN 

CPA/S_CDF3_NONLIN 

CP  A/S_PDF_GT_N  ONLI 

M1PA/SERIALNONLIN 

LINEAR/SERIALLIN 


Description 

Histogram  MNP 
EBHMNP 
FSA  MNP 

CPA  with  Linear  Transform  MNP 

CPA  with  Discontinuous  Auxiliary  Function  MNP 

CPA  with  Gaussian  Tails  MNP 

MIPA  MNP  (for  20d  and  4*  order  MIPA) 

Linear  receiver 


NONLINEAR  RECEIVERS  WITHOUT  CORRELATOR: 


Librarv/Filename 

fflSTO/S_NONLIN_NC 

HISTO/S_ENONLIN_NC 

FOURIER/SNONLINNC 

CPA/S_CDF2_NL_NC 

CPA/S_CDF3_NLNC 

CPA/S_PDF_GT_NL_NC 

MIPA/S_NONUN_NC 


PSSCrip.tiOQ 
Histogram  MNP 
EBHMNP 
FSA  MNP 

CPA  with  Linear  Transform  MNP 

CPA  with  Discontinuous  Auxiliary  Function  MNP 

CPA  with  Gaussian  Tails  MNP 

MIPA  MNP  (for  2nd  and  4*  order  MIPA) 


Table  1 


1  In  this  report  SPW  filenames  are  written  as  LIBRARY/FILENAME,  where  LIBRARY  is  the 
Block  Diagram  Editor  (BDE)  library  in  which  the  block  or  system  is  stored,  and 
FILENAME  is  the  name  of  the  block  or  system. 
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There  are  two  main  categories  of  MNPs:  MNPs  with  a  correlator  and  MNPs  without  a 
correlator.  The  MNPs  with  a  correlator  implement  the  entire  memoryless  Locally  Optimal  (LO) 
algorithm  [Illi93b]: 

Choose  (4,4,)  which  maximizes: 

K^cos^^gCr^sin^}  (1) 


where  g(rt)=- 


/*(r*) 


is  the  LO  Memory  less  Nonlinear  Transform  (MNT). 


Since  the  correlator  performs  a  decision,  these  blocks  constitute  the  entire  nonlinear  receiver, 
and  may  be  used  in  systems  such  as  the  conventional  Quadrature  Phase  Shift  Keying  (QPSK) 
[Taub86]  systems  to  be  discussed  in  Section  3.  It  is  important  to  make  the  distinction  between 
the  MNT,  which  is  the  mathematical  formula,  and  the  MNP,  which  is  the  block  that  applies  the 
MNT. 


The  MNPs  without  a  correlator  only  implement  the  MNT  and  convert  it  to  rectangular 
coordinates  as  follows: 

output  =  [g,t  sCi]r = [sL g(rj  cos0t  s^gir,)  sindf  (2) 

These  blocks  may  be  used  in  systems  where  additional  post-processing  is  required  after  the 
MNP,  as  in  the  Spread  Spectrum  (SS)  systems  to  be  discussed  in  Section  3. 
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3.  The  ITT  SPW  Systems 


Nine  modular  SPW  systems  have  been  designed  to  incorporate  the  MNPs  previously 
discussed.  Four  of  these  systems  are  conventional  QPSK  systems  and  five  are  Direct  Sequence 
(DS)  SS  [Taub86]  QPSK  systems.  The  conventional  systems  aid  in  isolating  and  characterizing 
the  performance  of  the  MNPs,  and  the  DSSS  systems  help  analyze  the  performance  of  the  MNPs 
in  a  spread  spectrum  environment.  Both  the  conventional  and  the  DSSS  systems  are  based  on 
QPSK  modulation;  however,  it  would  be  straightforward  to  modify  them  to  employ  alternative 
modulation  schemes  such  as  Minimum  Shift  Keying,  Quadrature  Amplitude  Modulation,  or  M- 
ary  Phase  Shift  Keying.  The  filenames  and  purposes  of  the  IIT  SPW  systems  are  summarized 
in  Table  2.  All  systems  simulate  a  communication  system  with  jammers  and  nonlinear 
processing,  and  each  one  has  certain  processing  after  the  nonlinearity  that  is  also  listed  in  the 
table. 


Figure  (1)  shows  the  basic  conventional  QPSK  system,  SIM/PE,  which  calculates  the 
probability  of  bit  error  (Pb)  for  any  of  the  serial  MNPs  or  the  linear  receiver.  Since  QPSK  is 
a  basic  form  of  modulation  it  was  chosen  to  isolate  the  effect  of  the  MNP  on  the  overall  system 
performance.  The  rl/complex_data  block  generates  random  QPSK  data,  and  the 
JAM/CHANNEL  block  adds  the  jammers  and  Gaussian  noise.  An  MNP  or  the  linear  receiver 
block  must  be  copied  by  the  user  to  the  designated  space  in  Fig.  (1).  The 
histo/prob_of_error  block  counts  the  number  of  errors  the  receiver  makes  in  decoding  the 
bits  and  computes  Pb. 

The  sim/PES  system  shown  in  Fig.  (2)  is  able  to  calculate  Pb  for  up  to  three  different 
receiver  configurations  simultaneously  to  facilitate  comparisons  of  the  MNPs.  The  SIM/PE  and 
SIM/PES  simulations  may  be  used  in  conjunction  with  the  SPW  iteration  macro  [Vam93]  to 
generate  Bit  Error  Rate  (BER)  curves  as  a  function  of  various  system  parameters. 
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Filenames  of  ITT  SPW  Systems 


CONVENTIONAL  QPSK  SIMULATIONS: 


Utaarxffikname 

SIM/PE 

SIM/PES 

SIM/HSIGS 

SIM/PSIGS 


DsSfiOBliQn 

Computes  Pb  for  any  MNP 

Computes  Pb  for  any  three  MNPs 

Computes  Pb  &  plots  signals  for  histogram  MNPs 

Computes  Pb  &  plots  signals  for  polynomial  and  FSA  MNPs 


QPSK  SPREAD  SPECTRUM  SIMULATIONS: 


Librarv/Filename 

SIM/SS_PE 

SIM/SSPES 

SIM/DSSSH 

SIM/DSSSP 

SIM/DSSSL 


Description 

Computes  Pb  for  any  MNP 

Computes  Pb  for  any  three  MNPs 

Computes  Pb  &  plots  signals  for  histogram  MNPs 

Computes  Pb  &  plots  signals  for  polynomial  and  FSA  MNPs 

Computes  Pb  &  plots  signals  for  linear  system 


Table  2 

The  other  two  conventional  systems  are  SIM/HSIGS  and  sim/psigs,  and  are  displayed  in 
Figs.  (3)  and  (4),  respectively.  HSIGS  is  an  abbreviation  of  "histogram  signals".  In  addition  to 
computing  Pb,  this  system  outputs  all  relevant  signals  for  detailed  analysis.  The  SIM/HSIGS 
system  works  for  both  the  histogram  and  EBH  MNPs.  PSIGS  is  short  for  "polynomial  signals", 
and  the  SIM/PSIGS  system  works  for  the  CPA  and  MIPA  MNPs.  sim/psigs  also  works  for  the 
FSA  MNP,  even  though  the  FSA  is  not  a  polynomial  approximation  method. 

The  DSSS  systems  are  shown  in  Figs.  (5)  through  (9).  The  basic  system,  sim/ss_pe,  is 
shown  in  Fig.  (5).  QPSK  data  is  generated  by  the  SPW  QPSK  SOURCE  library  block.  The 
spectral  spreading  is  achieved  by  multiplying  the  QPSK  data  by  a  Pseudo  Noise  (PN)  sequence, 
and  the  jammers  and  Gaussian  noise  are  added  in  the  channel.  The  user  must  copy  an  MNP 
without  a  correlator  into  the  designated  space.  The  transformed  sequence  is  multiplied  by  the 


5 


same  PN  sequence  to  despread  the  spectrum.  A  PSK  matched  filter  demodulator  library  block 
is  used  to  receive  the  message,  and  the  HISTO/PROB_OF_ERROR  block  is  used  to  compute  P». 


4.  Using  the  MNPs  in  an  UT  SPW  System 

A  simple  three  step  process  is  required  to  use  the  MNPs  in  the  IIT  SPW  simulations. 
First,  the  user  enters  the  Block  Diagram  Editor  (BDE)  and  loads  one  of  the  systems  listed  in 
Table  2  into  a  viewport.  Second,  the  desired  MNP  is  copied  from  the  User  Palette  into  the 
designated  space  in  the  system,  making  sure  that  all  wires  are  properly  connected.  The 
parameters  in  the  MNP  block  are  already  exported  to  the  parameters  on  the  top  level  of  the 
system.  Finally,  the  system  level  parameters  are  adjusted  by  the  user  to  their  desired  settings 
and  the  system  is  simulated. 


5.  Parameters  for  IIT  SPW  Simulations 

There  are  numerous  parameters  which  specify  the  configuration  of  the  IIT  SPW 
simulations.  It  is  not  possible  to  specify  every  parameter  explicitly  since  some  parameters  are 
dependent  on  others.  For  example,  the  symbol  rate,  R„  the  sampling  frequency,  f,  and  the 
number  of  samples  per  symbol,  N„  are  related  by  R,  -fJNr  Thus,  specifying  any  two  of  these 
parameters  implies  a  value  for  the  third.  In  this  manual  a  parameter  which  is  directly  specified 
is  enclosed  in  braces  {}  for  clarity.  The  remaining  parameters  are  considered  to  be  variables , 
and  their  values  are  a  function  of  the  specified  parameters.  A  listing  of  the  parameters  of  the 
IIT  SPW  simulations  is  presented  in  Table  3. 
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ITT  SPW  System  Parameters 


5.1  Parameter  Overview 


Several  of  the  parameters  in  the  SPW  simulations  are  interrelated,  and  certain 
relationships  between  them  must  be  maintained.  In  addition,  the  number  of  iterations  per 
simulation  run  is  also  dependent  on  certain  parameters.  It  is  necessary  to  choose  the  correct 
number  of  iterations  per  simulation  to  be  certain  that  the  desired  number  of  Pb  values  will  be 
calculated.  Do  not  select  run  to  EOF,  because  there  is  no  end  of  file  in  these  systems  and  they 
would  run  indefinitely. 

The  nonlinearity  in  the  MNP  blocks  requires  vector  operation.  In  all  of  the  MNPs  listed 
in  Table  1  the  incoming  serial  data  is  buffered,  operated  on  by  the  nonlinearity  in  vector  form, 
and  converted  back  to  serial  form.  The  length  of  this  vector  is  N,  the  number  of  samples  per 
correlation.  This  length  also  corresponds  to  the  number  of  samples  for  each  MNT 
approximation.  The  correlation  must  be  based  on  an  integer  number  of  symbols,  so  N/N,  must 
be  an  integer.  If  N/N,  is  not  an  integer  correlations  after  the  first  one  will  not  be  synchronized 
to  the  symbol  period  and  Pb  will  suffer  as  a  result. 

Let  IQ  be  the  number  iterations  for  a  conventional  QPSK  simulation  run  and  k  be  the 
number  of  Pb  points  to  be  computed.  Also,  let  C  be  the  number  of  correlations  per  Pb 
calculation.  This  also  corresponds  to  the  number  of  MNT  approximations  per  Pb  calculation. 
C  is  an  implicit  variable  which  is  a  function  of  other  parameters.  The  following  three 
constraints ,  (A),  (B),  and  (C),  must  be  maintained  at  all  times: 


N/N,*  C  =  Nn 

(A) 

lQ-k*N*C  +  N 

(B) 

N/N,  must  be  an  integer 

(Q 

Fundamental  Constraint  Box  (1) 
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For  example,  let  {N}  =  5,000,  {N,}  =  25  and  {Nn}  =  1,000.  (These  are  the  default  values 

in  the  IIT  SPW  systems.)  From  constraint  (A), 

{N)/{Nt}*C={Nn} 

5,000/25  *  C=  1,000 
C=5 

The  number  of  iterations  required  for  a  single  Pb  data  point  (k  =  1)  is  determined  from 
constraint  (B)  to  be  IQ  =  25,000  +  5,000  =  30,000.  To  compute  three  Pb  data  points  requires 
IQ  =  3*  25,000  +  5,000  =  80,000.  Note  that  for  both  cases  C  =  5  and  {/V}  /  {N,}  =  200, 
which  is  an  integer  as  required  by  constraint  (C).  When  using  the  SPW  iteration  macro  to 
iterate  over  N„  N  must  be  an  integer  multiple  of  the  least  common  multiple  of  all  the  values  of 
Nt.  This  insures  that  constraint  (C)  is  satisfied  throughout  the  iteration. 

The  additional  N  term  in  constraint  (B)  is  a  result  of  buffering  requirements.  If  the 
number  of  iterations  for  each  simulation  are  chosen  to  be  greater  than  the  required  number 
more  than  k  Pb  values  may  be  computed.  Conversely,  if  the  number  of  iterations  is  less  than 
required,  fewer  than  k  Pb  values  will  be  computed.  For  the  Spread  Spectrum  systems  in  this 
report,  the  required  number  of  iterations  is  Is  —  IQ  +  N„  due  to  time  delay  in  the  PSK 
demodulator  block. 

Note:  The  CPA  with  Gaussian  Tails 
nonlinearities  (CPA/PDF_GT_NONLI  and 
CPA/PDF_GT_NL_NC)  must  have  at  least  four 
bins  to  function  properly. 

Fundamental  Constraint  Box  (2) 
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5.2  Detailed  Parameter  Discussion 


This  following  discussion  outlines  the  reasons  why  the  parameters  in  Table  3  are  the 
preferred  way  to  specify  the  configuration  of  the  SPW  systems.  In  addition,  two  very  important 
relationships  which  are  summarized  below  are  discussed: 


Maintain  at  ail  times 

For  the  SS  systems,  also  maintain  N^NJKj  at  all  times 


Fundamental  Constraint  Box  (3) 

The  derivations  in  the  following  sections  are  performed  in  each  channel  separately. 
QPSK  modulation  may  be  viewed  as  Binary  Phase  Shift  Keying  (BPSK)  in  the  In-Phase  and 
Quadrature  channels;  each  QPSK  symbol  consists  of  two  BPSK  bits.  The  Continuous  Wave 
(CW)  jammer  may  be  viewed  as  a  sinusoid  in  each  channel,  and  the  Partial  Band  (PB)  jammer 
may  be  viewed  as  a  filtered  Gaussian  signal  in  each  channel.  EJN0  and  J/S  are  the  same  for 
each  channel  and  for  the  composite  signal.  In  addition,  the  bit  rate  of  each  channel,  /^,  is  equal 
to  R„  and  the  number  of  samples  per  bit,  Nb,  is  equal  to  Ns. 

5.2.1  System  Parameters  Related  to  the  Information  Signal  and  Gaussian  Noise 
Signal  Energy 

The  energy,  E,  of  a  signal  s(t)  is  defined  as 

m 

EaJ|s(r)|2dr  (3) 

-OB 

For  discrete  signals  this  becomes 
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(4) 


E-t  |J(")I!7 

Jg 


where  f  is  the  sampling  frequency.  A  BPSK  bit  maintains  a  constant  level  ±A  for  the  entire 
bit  period.  The  energy  of  one  bit,  Ebt  is  given  by 


E-V  —  =  -  — 

b~h  f,a  f,  '  r. 


f 

where  Ns  is  equal  to  the  number  of  samples  per  bit  and  Rt*—  is  equal  to  the  bit  rate. 

Nt 


(5) 


In  many  analyses  it  is  desirable  to  compute  Pb  as  a  function  of  EJN0,  where  N0  is  the 
background  noise  power.  This  can  be  achieved  by  making  Eb/N0  a  parameter,  {E^/A/q}  ,  and 
making  N0  a  parameter,  {W0} .  Thus,  for  BPSK 


(w„)  (W0}R, 


(6) 


Then  the  bit  amplitude  is 


A2 


=  {VA« 


(7) 


gjgnql  PQwgr 

The  power,  P,  of  a  signal  is  defined  as 

where  T  is  the  period  of  the  signal.  For  discrete  signals  the  power  is 
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(9) 


w«)i*7 

1  T  Jt 

For  a  BPSK  signal  with  amplitude  ±A  the  power,  5,  is 

— -/l2 

"s&f. 


5 2,1  System  Parameters  Relating  the  Continuous  Wave  Jammer  to  the  Information  Signal 


The  power  of  a  CW  jammer  is 


2”  j AIc°s2(<aJt)dt * y- 


where  Aj  is  the  jammer  amplitude  and  <u;  is  the  jammer  frequency.  The  jammer  amplitude  can 
be  found  as  a  function  of  the  parameter  {.//#„}  : 

um.j-.da  a® 


A?  *2{IIN0)(N0) 


It  is  often  more  useful  to  write  the  jammer  amplitude  as  a  function  of  the  parameter  {//s} . 

UlA.lmdl  a- 

S  A2 
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(15) 


/t/.24!M  -2 {£,«„} i^l//s} 


Relationship  between  CW  Jammer  Parameters 
It  is  quite  straightforward  to  show  that 


{/«} 


(£*/w4 


(1«) 


Usually  these  parameters  are  specified  in  deciBels  (dB),  and  the  following  relationship  results: 


N. 


(17) 


Note  that  in  Eq.  (17)  [EbIN0}  is  not  in  dB. 


Choice  of  CW  Jammer  Parameters 

One  of  the  key  assumptions  in  the  LO  derivation  is  that  the  signal  is  small  compared  to 
the  jammer,  i.e.  the  ratio  A,2 /A2  is  large.  From  Eqs.  (7)  and  (13)  it  is  seen  that  if  the  parameter 
{J/N0}  is  used,  this  ratio  becomes 


a}  _  2{ffly„} 

s 


(18) 


This  means  that  the  validity  of  the  small  signal  assumption  is  not  only  dependent  on  {J/N0},  but 
also  on  many  other  parameters.  However,  if  the  parameter  {J/S}  is  used  instead  of  {J/N0},  the 
ratio  A2 /A2  is  determined  from  Eqs.  (7)  and  (15)  to  be 
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(1») 


A>  2{£,/N0}i^lll//, 


fr'No) 


Now  the  small  signal  assumption  is  only  dependent  on  {//£}.  For  this  reason,  the  use  of  the 
{7/5}  parameter  is  preferred  over  {J/N0} . 


The  message  bandwidth  is  proportional  to  the  bit  rate,  and  R _»-i.  If  N,  is  varied,  the 
message  bandwidth  will  change  but  fJt  the  jammer  frequency,  will  not.  Thus,  the  spectrum  of 
the  message  sequence  will  expand  or  contract  while  the  spectral  position  of  the  jammer  will 
remain  unchanged.  As  the  relative  position  of  the  jammer  and  the  message  change,  the 
performance  of  the  system  will  change  dramatically.  In  order  to  isolate  the  performance  of  the 
MNP  from  the  effects  of  this  shift,  it  is  useful  to  use  a  frequency  ratio,  {£//?,},  instead  of  {/)}. 
Using  this  ratio,  the  jammer  frequency  is  f-{fJR,}R,-  This  way  fj  and  the  message  bandwidth 
change  in  the  same  manner  when  Rt  is  varied. 

Alternatively,  the  sampling  frequency  may  be  set  to  f,={Ns}{R,}.  Then  if  R,  is  constant, 
the  message  bandwidth  will  remain  constant  as  N,  is  varied.  The  spectral  positions  of  both  the 
jammer  and  the  message  signal  will  remain  constant. 
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5.2.3  System  Parameters  Relating  the  Partial  Band  Jammer  to  the  Information  Signal 
PB  Jammer  Power 

The  autocorrelation  function,  rw(m),  for  sampled  white  Gaussian  noise,  vv(/i),  is  given  by 

rw(m)=aj{m)  <2°) 

where  a*  is  the  variance  of  w(n)  and  6(m)  =  |q’  is  the  Kronecker  delta  function.  The 
Power  Spectral  Density  (PSD),  Sw(f),  is  given  by 

S„(/Wi  foraU/  01) 


The  frequency  response  of  an  ideal  lowpass  filter  is 


Wf) 


lA<fc 

fc<\f\<\ 


(22) 


where  f.  is  the  cutoff  frequency  (normalized  by  the  sampling  frequency,^).  The  Gaussian  noise 
PB  jammer  is  constructed  by  passing  white  Gaussian  noise  through  the  lowpass  filter  (for 
simulation  at  baseband).  Therefore,  the  PSD  of  the  PB  jammer,  SPB(f),  is 


1/1  </c 
fc<\f\<\ 


(23) 


The  power  of  the  PB  jammer,  JPB,  may  be  found  by  integrating  the  PSD  over  the  domain  of  /, 
resulting  in 
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{■ 

-1/2 


S„('a-f)df 


'V,°l 


Given  the  jammer-to-signal  ratio  parameter  [//*/$},  then 


2fA°l 

{EM  {*„}/, 


To  find  the  required  value  of  for  a  given  {7^/5} 


.2 .  tWA  {£*TO  TO  /, 

**'"  2  f.  N. 


{JPB/S}  {Eb/N0}  {N0} 

' — iz — • 


16 


If  fc  is  chosen  as  a  fraction  of  the  bit  rate,  i.e.  if  {£/#,}  is  used  then 


.  (Vi!  !W  R! 

"■  WM — 


07) 


5.2.4  Spread  Spectrum  System  Parameters 


The  processing  gain,  number  of  samples  per  chip,  and  number  of  samples  per  symbol 
are  related  by 


Ns  =  NcPG 


(28) 


Also,  since  R=^L,  then  the  chip  rate,  Re,  is  equal  to 

Ns 

f  N 

R  =LL=R  —=R  PG 
e  Nc  3  Nc  3 


(29) 
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Appendix  A 


Docgen  Listing 


This  Appendix  is  an  alphabetical  listing  of  the  help  screens  available  for  each  in  SPW 
block.  In  this  Appendix  the  inputs,  parameters  and  outputs  of  each  block  are  capitalized  in  the 
description  of  the  block  to  distinguish  them  from  the  text.  However,  in  the  listing  of  the  inputs, 
parameters  and  outputs,  the  capitalization  matches  the  capitalization  of  the  actual  parameter 
names  in  SPW. 

The  names  of  hierarchical  blocks  are  marked  with  a  dagger f.  The  details  of  these  blocks 
are  printed  in  alphabetical  order  in  Appendix  B. 


A-l 


Name: 

cpa/cdf2+ 

Description: 

This  hierarchical  block  generates  a  Cumulative  Distribution  Function  (CDF)  and  a  Probability 
Density  Function  (PDF)  of  the  magnitude  of  the  input  message  by  applying  the  Continuous 
Polynomial  Approximation  (CPA)  algorithm  to  a  histogram  CDF.  A  vec/heap_sort  block  is  used 
to  sort  the  data  to  prepare  it  for  the  histo/equi  block,  which  generates  the  histogram  CDF.  The 
cpa/slopel  block  computes  the  derivatives  of  the  histogram  needed  by  the  CPA  algorithm.  The 
cpa/coefjt  block  uses  the  equiprobable  bin  histogram  CDF,  its  derivatives,  and  the  breakpoints 
to  compute  the  CPA  CDF  and  PDF. 

This  block  applies  the  linear  transform  to  the  polynomial.  Without  the  linear  transform  the 
magnitude  values  in  each  bin  would  range  from  BP[K-1]  to  BP[K].  The  linear  transform  shifts 
each  bin  to  the  origin,  so  the  new  magnitude  values  range  from  0  to  BP[K]-BP[K-1]. 


Inputs: 

data 


Magnitude  of  input  message  sequence 


Parameters: 

bins  Number  of  bins 

samples  Number  of  samples 


Outputs: 

bp  Breakpoints  (CDF  and  PDF  interval  boundaries) 

cdf  CDF  polynomial  coefficients 

pdf  PDF  polynomial  coefficients 

See  also: 

vec/heap_sort,  histo/equi,  cpa/slopel,  cpa/coefjt 


Name: 

cpa/cdf2_nonlinf 

Description: 

This  is  a  hierarchical  block  which  applies  the  nonlinearity  to  the  input  signal  in  vector  format. 
The  histo/rec_to_polar  block  converts  the  two  dimensional  input  signal  from  rectangular  to  polar 
coordinates,  and  a  CPA  of  the  Probability  Density  Function  (PDF)  of  the  magnitude  is  computed 
by  the  cpa/cdf2  block.  The  poly/mnt  block  obtains  the  Memoryless  Nonlinear  Transform 


A-2 


(MNT)  from  the  CPA  PDF  and  applies  it  to  the  magnitude  of  the  input  message  sequence.  The 
result  is  converted  back  to  rectangular  coordinates  by  the  histo/polar_to_rec  block  using  the 
unmodified  phase.  The  poly/plot  block  generates  a  plot  of  the  PDF  and  CDF  from  their 
polynomial  coefficients. 

The  coefficients  in  the  cpa/cdf2  block  are  computed  with  the  linear  transform  applied,  thus 
each  bin  is  shifted  to  the  origin. 

Refer  to  the  poly/mnt  block  for  more  information  on  the  MNT. 


Inputs: 

i_in 

q_in 

In-Phase  component  of  input  message  sequence 
Quadrature  component  of  input  message  sequence 

Parameters: 

samples 

bins 

points 

Number  of  samples  in  input  vectors 

Number  of  bins  in  the  histogram  PDF 

Number  of  plot  points  for  PDF  and  MNT 

Outputs: 

pdf 

i_out 

q_out 

mnt 

bp 

CPA  PDF 

In-Phase  component  of  transformed  sequence 
Quadrature  component  of  transformed  sequence 
Plot  of  the  Memoryless  Nonlinear  Transform 
Breakpoints  (PDF  and  MNT  interval  boundaries) 

See  also: 

histo/rec_to_polar,  cpa/cdf2,  poly/mnt,  poly/plot,  histo/polar_to_rec 


Name: 

cpa/cdf3+ 

Description: 

This  hierarchical  block  generates  a  Cumulative  Distribution  Function  (CDF)  and  a  Probability 
Density  Function  (PDF)  of  the  magnitude  of  the  input  message  by  applying  the  Continuous 
Polynomial  Approximation  (CPA)  algorithm  to  a  histogram  CDF.  A  vec/heap_sort  block  is  used 
to  sort  the  data  to  prepare  it  for  the  histo/equi  block,  which  generates  the  histogram  CDF.  The 
cpa/slopel  block  computes  the  derivatives  of  the  histogram  needed  by  the  CPA  algorithm.  The 
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cpa/coefjtdaf  block  uses  the  equiprobable  bin  histogram  CDF,  its  derivatives,  and  the 
breakpoints  to  compute  the  CPA  CDF  and  PDF. 

This  block  (cpa/cdf3)  is  identical  to  the  cpa/cdf2  block  except  that  the  constraint  that  the  CDF 
be  continuous  has  been  relaxed.  Thus,  the  equations  to  determine  the  coefficient  values  are 
different. 

Like  the  cpa/cdf2  block,  this  block  applies  the  linear  transform  to  the  polynomial.  Without 
the  linear  transform  the  magnitude  values  in  each  bin  would  range  from  BP[K-1]  to  BP[K].  The 
linear  transform  shifts  each  bin  to  the  origin,  so  the  new  magnitude  values  range  from  0  to 
BP[K]-BP[K-1]. 


Inputs: 

data  Magnitude  of  input  message  sequence 


Parameters: 

bins  Number  of  bins 

samples  Number  of  samples 

Outputs: 

bp  Breakpoints  (CDF  and  PDF  interval  boundaries) 

cdf  CDF  polynomial  coefficients 

pdf  PDF  polynomial  coefficients 

See  also: 

vec/heap_sort,  histo/equi,  cpa/slopel,  cpa/coefjtdaf 


Name: 

cpa/cdf3_nonlin+ 

Description: 

This  is  a  hierarchical  block  which  applies  the  nonlinearity  to  the  input  signal  in  vector  format. 
The  histo/rec_to  jjolar  block  converts  the  two  dimensional  input  signal  from  rectangular  to  polar 
coordinates,  and  a  CPA  of  the  Probability  Density  Function  (PDF)  of  the  magnitude  is  computed 
by  the  cpa/cdf3  block.  The  poly/mnt  block  obtains  the  Memoryless  Nonlinear  Transform 
(MNT)  from  the  CPA  PDF  and  applies  it  to  the  magnitude  of  the  input  message  sequence.  The 
result  is  converted  back  to  rectangular  coordinates  by  the  histo/polar_to_rec  block  using  the 


A-4 


unmodified  phase.  The  poly/plot  block  generates  a  plot  of  the  PDF  and  CDF  from  their 
polynomial  coefficients. 

The  coefficients  in  the  cpa/cdO  block  are  computed  with  the  linear  transform  applied,  thus 
each  bin  is  shifted  to  the  origin.  Also,  the  constraint  that  the  CDF  be  continuous  is  removed. 

Refer  to  the  poly/mnt  block  for  more  information  on  the  MNT. 


Inputs: 

i_in 

q_in 

In-Phase  component  of  input  message  sequence 
Quadrature  component  of  input  message  sequence 

Parameters: 

samples 

bins 

points 

Number  of  samples  in  input  vectors 

Number  of  bins  in  the  histogram  PDF 

Number  of  plot  points  for  PDF  and  MNT 

Outputs: 

pdf 

i_out 

q_out 

mnt 

bp 

CPA  PDF 

In-Phase  component  of  transformed  sequence 
Quadrature  component  of  transformed  sequence 
Plot  of  the  Memoryless  Nonlinear  Transform 
Breakpoints  (PDF  and  MNT  interval  boundaries) 

See  also: 

histo/rec_to_polar,  cpa/cdO,  poly/mnt,  poly/plot,  histo/polar_to_rec 


Name: 

cpa/coef_gt 

Description: 

This  block  computes  the  Continuous  Polynomial  Approximation  (CPA)to  a  function.  In  this 
implementation,  the  CPA  is  an  approximation  of  a  Probability  Density  Function  (PDF). 

The  coefficients  in  this  routine  are  computed  with  the  linear  transform  applied,  thus  each  bin 
is  shifted  to  the  origin.  The  first  and  last  bins  are  a  Gaussian  tail  instead  of  a  polynomial.  The 
mean  and  variance  for  the  Gaussian  are  taken  to  be  the  mean  and  variance  of  the  received  signal. 

Note  that  the  value  for  BINS  must  be  set  to  a  value  greater  than  or  equal  to  4. 


A-5 


WARNING: 

This  program  modifies  the  histogram  PDF  input.  If  another  block  is  going  to  also  use  the 
histogram  PDF  results  will  be  unpredictable. 

Inputs: 
bp 
f 

df 

mean 
var 

Parameters: 

bins  Number  of  bins 

Outputs: 

pdf  CPA  PDF 


Name: 

cpa/coefjt 

Description: 

This  block  computes  the  Continuous  Polynomial  Approximation  (CPA)to  a  function.  In  this 
implementation,  the  CPA  is  an  approximation  of  a  Cumulative  Distribution  Function  (CDF). 

The  coefficients  in  this  routine  are  computed  with  the  linear  transform  applied,  thus  each  bin 
is  shifted  to  the. origin. 

Inputs: 

bp  Breakpoints  (CDF  and  PDF  interval  boundaries) 

histo_cdf  Histogram  CDF 

df  First  derivative  of  Histogram  CDF 

ddf  Second  derivative  of  Histogram  CDF 

Parameters: 

bins  Number  of  bins 

Outputs: 

cdf  CPA  CDF 

pdf  CPA  PDF 


Breakpoints  (PDF  interval  boundaries) 
Histogram  PDF 

First  derivative  of  Histogram  PDF 
Average  of  the  received  signal 
Variance  of  the  received  signal 
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Name: 

cpa/coefjtdaf 

Description: 

This  block  computes  the  Continuous  Polynomial  Approximation  (CPA)to  a  function.  In  this 
implementation,  the  CPA  is  an  approximation  of  a  Cumulative  Distribution  Function  (CDF). 


The  coefficients  in  this  routine  are  computed  with  the  linear  transform  applied,  thus  each  bin 
is  shifted  to  the  origin.  Also,  the  constraint  that  the  CDF  be  continuous  is  removed. 


Inputs: 

bp  Breakpoints  (CDF  and  PDF  interval  boundaries) 

df  First  derivative  of  Histogram  CDF 

ddf  Second  derivative  of  Histogram  CDF 


Parameters: 

bins 


Number  of  bins 


Outputs: 

cdf  CPA  CDF 

pdf  CPA  PDF 


Name: 

cpa/pdf_gtf 

Description: 

This  hierarchical  block  generates  a  Probability  Density  Function  (PDF)  of  the  magnitude  of 
the  input  message  by  applying  the  Continuous  Polynomial  Approximation  (CPA)  algorithm  to 
a  histogram  PDF.  A  vec/heap_sort  block  is  used  to  sort  the  data  to  prepare  it  for  the  histo/equi 
block,  which  generates  the  histogram  PDF.  The  vec/ave_bp  obtains  values  of  the  histogram  at 
the  breakpoints  by  averaging  the  values  of  the  adjacent  bins.  The  cpa/slope2  block  computes 
the  derivatives  of  the  histogram  needed  by  the  CPA  algorithm.  The  cpa/coef_gt  block  uses  the 
equiprobable  bin  histogram,  its  derivatives,  the  breakpoints,  and  the  mean  and  variance  of  the 
magnitude  to  compute  the  CPA  PDF. 

The  first  and  last  bins  are  approximated  using  Gaussian  tails  instead  of  polynomials.  Because 
of  this,  the  poly/mnt _gt  and  poly/plot_gt  blocks  should  be  used  with  this  block. 
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'Hie  linear  transform  is  applied  to  the  polynomial  bins.  Without  the  linear  transform  the 
magnitude  values  in  each  bin  would  range  from  BP[K-1]  to  BP[K].  The  linear  transform  shifts 
each  bin  to  the  origin,  so  the  new  magnitude  values  range  from  0  to  BP[K]-BP[K-1]. 

Note  that  the  value  for  BINS  must  be  set  to  a  value  greater  than  or  equal  to  4. 


Inputs: 

data 


Magnitude  of  input  message  sequence 


Parameters: 

bins  Number  of  bins 

samples  Number  of  samples 

Outputs: 

bp  Breakpoints  (PDF  interval  boundaries) 

pdf  PDF  polynomial  coefficients 

See  also: 

vec/heap_sort,  histo/equi,  cpa/slope2,  cpa/coefjtdaf,  poly/mnt_gt,  poly/plot_gt 


Name: 

cpa/pdf_gt_nonlinf 

Description: 

This  is  a  hierarchical  block  which  applies  the  nonlinearity  to  the  input  signal  in  vector  format. 
The  histo/rec_to_polar  block  converts  the  two  dimensional  input  signal  from  rectangular  to  polar 
coordinates,  and  a  CPA  of  the  Probability  Density  Function  (PDF)  of  the  magnitude  is  computed 
by  the  cpa/pdf_gt  block.  The  poly/mnt  block  obtains  the  Memoryless  Nonlinear  Transform 
(MNT)  from  the  CPA  PDF  and  applies  it  to  the  magnitude  of  the  input  message  sequence.  The 
result  is  converted  back  to  rectangular  coordinates  by  the  histo/polar_to_rec  block  using  the 
unmodified  phase.  The  poly/plot  block  generates  a  plot  of  the  PDF  and  CDF  from  their 
polynomial  coefficients. 

The  coefficients  in  the  cpa/pdf_gt  block  are  computed  with  the  linear  transform  applied,  thus 
each  bin  is  shifted  to  the  origin.  The  first  and  last  bins  are  a  Gaussian  tail  instead  of  a 
polynomial.  The  mean  and  variance  for  the  Gaussian  are  taken  to  be  the  mean  and  variance  of 
the  received  signal. 

Refer  to  the  poly/mnt  block  for  more  information  on  the  MNT. 
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Inputs: 

ijn 

In-Phase  component  of  input  message  sequence 

q_in 

Quadrature  component  of  input  message  sequence 

Parameters: 

samples 

Number  of  samples  in  input  vectors 

bins 

Number  of  bins  in  the  histogram  PDF 

points 

Number  of  plot  points  for  PDF  and  MNT 

Outputs: 

pdf 

CPA  PDF 

i_out 

In-Phase  component  of  transformed  sequence 

q_°ut 

Quadrature  component  of  transformed  sequence 

mnt 

Plot  of  the  Memoryless  Nonlinear  Transform 

bp 

Breakpoints  (PDF  and  MNT  interval  boundaries) 

See  also: 

histo/rec_to_polar. 

cpa/pdf_gt,  poly/mnt,  poly/plot,  histo/polar_to_rec 

Name: 

cpa/s_cdf2_nl_ncf 

Description: 

This  hierarchical  block  buffers  the  received  data  and  converts  it  from  serial  to  vector  data  and 
applies  it  to  the  cpa/cdf2_nonlin  nonlinearity.  This  block  does  not  generate  the  summation  part 
of  the  LO  decision  statistic. 

The  romelib/timing  block  generates  timing  waveforms  which  enable  and  disable  the 
nonlinearity  and  all  subsequent  blocks  at  appropriate  times.  Of  particular  interest  are  the 
CLK_OUT  and  HOLD_VEC  signals.  The  HOLD_VEC  signal  holds  vector  blocks  while  the 
serial  data  is  buffered  into  vector  form.  This  signal  goes  low  during  the  simulation  iteration  in 
which  the  vector  processing  completes  and  the  outputs  are  available.  The  CLK_OUT  signal 
provides  flexibility  for  use  in  multirate  systems  by  pulsing  high  when  a  valid  data  sample  is 
produced  at  the  output  of  the  cpa/s_cdf2_nl_nc  block.  The  CLK_OUT  signal  has  the  same 
period  as  the  CLK_IN  signal,  thus  output  samples  are  produced  at  the  same  rate  as  that  of  the 
input  signal.  The  remaining  romelib/timing  signals  are  internal  to  the  cpa/s_cdf2_nl_nc  block, 
and  are  used  in  serial  to  vector  and  vector  to  serial  buffering. 
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The  SAMPLES_PER_SYMBOL,  NUM_REFS,  and  M_TYPE  are  characteristics  of  the 
reference  symbols.  The  reference  signals  must  be  created  before  using  this  block,  and  are  stored 
in  the  /spwdata/ref  library.  This  permits  the  designer  to  use  any  desired  signal  constellation. 
The  NUM_REFS  parameter  must  be  exported  to  the  histo/correlator2  block,  and  the  other  two 
parameters  are  exported  to  the  Filename  of  the  rl/inf_vsource  block  as  "ref/"&substr 
("m_type":model,2,length("m_type":model)-2)&"_"&xstring(Msamples_per_symbor:model)& 
"spb_i,H  for  the  I  channel  signals  and  " ref/ "  &substr( " m_type " : model , 2*, length( ” m_type " : 
model)~2)  &"_"&xstring("samples_per_symbor‘:model)&Mspb_q”'  for  the  Q  channel  signals. 
For  example,  the  filenames  for  a  Binary  Phase  Shift  Keying  (BPSK)  signal  with  25  samples  per 
symbol  would  be  /spwdata/ref/bpsk_25spb_i  and  /spwdata/ref/bpsk_25spb_q. 


Inputs: 

in 

Input  message  sequence  (complex  input) 

clk_in 

Clock  input  for  romeiib/timing 

Parameters: 

samples 

Number  of  samples  in  input  vectors 

bins 

Number  of  bins  in  the  Probability  Density  Function  (PDF) 

samples_per_symbol 

Number  of  samples  per  data  symbol 

num_refs 

Number  of  reference  signals 

m_type 

Modulation  type 

Outputs: 

mnt 

Plot  of  the  Memoryless  Nonlinear  Transform 

g 

Complex  transformed  sequence  (vector) 

out 

Complex  transformed  sequence  (serial) 

clk_out 

Timing  signal:  goes  high  when  valid  data  sample  is  available  at  the  output 
of  cpa/s_cdf2_nl_nc 

hold^vec 

Timing  signal:  goes  low  when  vector  outputs  are  available 

pdf 

Plot  of  the  PDF  of  magnitude 

bp 

Breakpoints  (PDF  and  MNT  interval  boundaries) 

See  also: 

cpa/cdf2_nonlin, 

romeiib/timing,  rl/inf_vsource,  cpa/s_cdf2_nonlin 
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Name: 

cpa/s_cdf2_nonlinf 

Description: 

This  hierarchical  block  buffers  the  received  data  and  converts  it  from  serial  to  vector  data  and 
applies  it  to  the  cpa/cdf2_nonlin  nonlinearity.  In  addition,  the  summation  in  the  LO  decision 
statistic  is  formed  using  the  histo/correlator2  block.  The  output  of  the  correlator  is  converted 
back  to  serial  form. 

The  romelib/timing  block  generates  timing  waveforms  which  enable  and  disable  the 
nonlinearity  and  all  subsequent  blocks  at  appropriate  times.  Of  particular  interest  are  the 
CLK_OUT  and  HOLD_VEC  signals.  The  HOLD_VEC  signal  holds  vector  blocks  while  the 
serial  data  is  buffered  into  vector  form.  This  signal  goes  low  during  the  simulation  iteration  in 
which  the  vector  processing  completes  and  the  outputs  are  available.  The  CLK_OUT  signal 
provides  flexibility  for  use  in  multirate  systems  by  pulsing  high  when  a  valid  data  sample  is 
produced  at  the  output  of  the  cpa/s_cdf2_nonlin  block.  The  CLK_OUT  signal  has  the  same 
period  as  the  CLK_IN  signal,  thus  output  samples  are  produced  at  the  same  rate  as  that  of  the 
input  signal.  The  remaining  romelib/timing  signals  are  internal  to  the  cpa/s_cdf2_nonlin  block, 
and  are  used  in  serial  to  vector  and  vector  to  serial  buffering. 

The  SAMPLES_PER_SYMBOL,  NUM_REFS,  and  M_TYPE  are  characteristics  of  the 
reference  symbols.  The  reference  signals  must  be  created  before  using  this  block,  and  are  stored 
in  the  /spwdata/ref  library.  This  permits  the  designer  to  use  any  desired  signal  constellation. 
The  NUM_REFS  parameter  must  be  exported  to  the  histo/corre!ator2  block,  and  the  other  two 
parameters  are  exported  to  the  filename  of  the  rl/inf_vsource  block  as  "ref/"&substr 

("m_type":m°del,2,length("m_type":model)-2)&H_"&xstring(,,samples_per_symbor:model)& 

"spb_i’M  for  the  I  channel  signals  and  Href/"&substr(Hm_type":model,2,length(Hm_type": 
model)-2)  &"_"&xstring("samples_per_symbor,:model)&"spb_q’n  for  the  Q  channel  signals. 
For  example,  the  filenames  for  a  Binary  Phase  Shift  Keying  (BPSK)  signal  with  25  samples  per 
symbol  would  be  /spwdata/ref/bpsk_25spb_i  and  /spwdata/ref/bpsk_25spb_q. 

Inputs: 

in  Input  message  sequence  (complex  input) 

clk_in  Clock  input  for  romelib/timing 

Parameters: 

samples  Number  of  samples  in  input  vectors 

bins  Number  of  bins  in  the  Probability  Density  Function  (PDF) 

samples_per_symbol  Number  of  samples  per  data  symbol 

num_refs  Number  of  reference  signals 

m_type  Modulation  type 
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Outputs: 

mnt  Plot  of  the  Memoryless  Nonlinear  Transform 

g  Complex  transformed  sequence 

out  Decision  of  nonlinearity 

clk_out  Timing  signal:  goes  high  when  valid  data  sample  is  available  at  the  output 

of  cpa/s_cdf2_noniin 

hold_vec  Timing  signal:  goes  low  when  vector  outputs  are  available 

pdf  Plot  of  the  PDF  of  magnitude 

bp  Breakpoints  (PDF  and  MNT  interval  boundaries) 


See  also: 

cpa/cdf2_nonlin,  histo/correlator2,  romelib/timing,  rl/inf_vsource 


Name: 

cpa/s_cdf3_nl_ncf 

Description: 

This  hierarchical  block  buffers  the  received  data  and  converts  it  from  serial  to  vector  data  and 
applies  it  to  the  cpa/cdf3_nonlin  nonlinearity.  This  block  does  not  generate  the  summation  part 
of  the  LO  decision  statistic. 

The  romelib/timing  block  generates  timing  waveforms  which  enable  and  disable  the 
nonlinearity  and  all  subsequent  blocks  at  appropriate  times.  Of  particular  interest  are  the 
CLK_OUT  and  HOLD_VEC  signals.  The  HOLD_VEC  signal  holds  vector  blocks  while  the 
serial  data  is  buffered  into  vector  form.  This  signal  goes  low  during  the  simulation  iteration  in 
which  the  vector  processing  completes  and  the  outputs  are  available.  The  CLK_OUT  signal 
provides  flexibility  for  use  in  multirate  systems  by  pulsing  high  when  a  valid  data  sample  is 
produced  at  the  output  of  the  cpa/s_cdf3_nl_nc  block.  The  CLK_OUT  signal  has  the  same 
period  as  the  CLK_IN  signal,  thus  output  samples  are  produced  at  the  same  rate  as  that  of  the 
input  signal.  The  remaining  romelib/timing  signals  are  internal  to  the  cpa/s_cdf3_nl_nc  block, 
and  are  used  in  serial  to  vector  and  vector  to  serial  buffering. 

The  SAMPLES_PER_SYMBOL,  NUM_REFS,  and  M_TYPE  are  characteristics  of  the 
reference  symbols.  The  reference  signals  must  be  created  before  using  this  block,  and  are  stored 
in  the  /spwdata/ref  library.  This  permits  the  designer  to  use  any  desired  signal  constellation. 
The  NUM_REFS  parameter  must  be  exported  to  the  histo/correlator2  block,  and  the  other  two 
parameters  are  exported  to  the  filename  of  the  rI/inf_vsource  block  as  "ref/"&substr 
("  m_type" :  model  ,2,length("  mjype"  :model)-2)&M_"  &xstring("  samples_per_symbol"  :model)& 


"spb_i’"  for  the  I  channel  signals  and  "ref/"&substr("m_type":model,2,length(“m_type": 
model)-2)  & M  _M  &xstring("  samples jper_symbol  " :  model)& "  spb_q’ "  for  the  Q  channel  signals. 
For  example,  the  filenames  for  a  Binary  Phase  Shift  Keying  (BPSK)  signal  with  25  samples  per 
symbol  would  be  /spwdata/ref/bpsk_25spb_i  and  /spwdata/ref/bpsk_25spb_q. 

Inputs: 

in  Input  message  sequence  (complex  input) 

clkjn  Clock  input  for  romelib/timing 

Parameters: 
samples 
bins 

samples_per_symbol 

numjrefs 
m_type 

Outputs: 

mnt  Plot  of  the  Memoryless  Nonlinear  Transform 

g  Complex  transformed  sequence  (vector) 

out  Complex  transformed  sequence  (serial) 

clk_out  Timing  signal:  goes  high  when  valid  data  sample  is  available  at  the  output 

of  cpa/s_cdf3_nl_nc 

hold  vec  Timing  signal:  goes  low  when  vector  outputs  are  available 

pdf  Plot  of  the  PDF  of  magnitude 

bp  Breakpoints  (PDF  and  MNT  interval  boundaries) 


See  also: 

cpa/cdf3_nonlin,  romelib/timing,  rl/inf_vsource,  cpa/s_cdO_nonlin 


Number  of  samples  in  input  vectors 

Number  of  bins  in  the  Probability  Density  Function  (PDF) 

Number  of  samples  per  data  symbol 

Number  of  reference  signals 

Modulation  type 
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Name: 

cpa/s_cdf3_nonlinf 

Description: 

This  hierarchical  block  buffers  the  received  data  and  converts  it  from  serial  to  vector  data  and 
applies  it  to  the  cpa/cdf3_nonlin  nonlinearity.  In  addition,  the  summation  in  the  LO  decision 
statistic  is  formed  using  the  histo/correlator2  block.  The  output  of  the  correlator  is  converted 
back  to  serial  form. 

The  Fomelib/timing  block  generates  timing  waveforms  which  enable  and  disable  the 
nonlinearity  and  all  subsequent  blocks  at  appropriate  times.  Of  particular  interest  are  the 
CLK_OUT  and  HOLD_VEC  signals.  The  HOLD_VEC  signal  holds  vector  blocks  while  the 
serial  data  is  buffered  into  vector  form.  This  signal  goes  low  during  the  simulation  iteration  in 
which  the  vector  processing  completes  and  the  outputs  are  available.  The  CLK_OUT  signal 
provides  flexibility  for  use  in  multirate  systems  by  pulsing  high  when  a  valid  data  sample  is 
produced  at  the  output  of  the  cpa/s_cdf3_nonlin  block.  The  CLK  OUT  signal  has  the  same 
period  as  the  CLK_IN  signal,  thus  output  samples  are  produced  at  the  same  rate  as  that  of  the 
input  signal.  The  remaining  romelib/timing  signals  are  internal  to  the  cpa/s_cdf3_nonlin  block, 
and  are  used  in  serial  to  vector  and  vector  to  serial  buffering. 

The  SAMPLES_PER_SYMBOL,  NUM_REFS,  and  M_TYPE  are  characteristics  of  the 
reference  symbols.  The  reference  signals  must  be  created  before  using  this  block,  and  are  stored 
in  the  /spwdata/ref  library.  This  permits  the  designer  to  use  any  desired  signal  constellation. 
The  NUMJREFS  parameter  must  be  exported  to  the  histo/correlator2  block,  and  the  other  two 
parameters  are  exported  to  the  filename  of  the  rl/inf_vsource  block  as  "ref/"&substr 
( "  m_type" :  model  ,2,length("  m_type” :  model)-2)&"_"  &xstring("samples _per_symbol" :  model)& 
"spb_i’n  for  the  I  channel  signals  and  Href/"&substr("m_type":model,2,length("m_type": 
model)-2)  &"_,,&xstring("samplesJ5er_symbo^:model)&"spb_q,"  for  the  Q  channel  signals. 
For  example,  the  filenames  for  a  Binary  Phase  Shift  Keying  (BPSK)  signal  with  25  samples  per 
symbol  would  be  /spwdata/ref/bpsk_25spb_i  and  / sp wdata/ ref/bpsk_25 spb_q . 

Inputs: 

in  Input  message  sequence  (complex  input) 

clk_in  Clock  input  for  romelib/timing 

Parameters: 

samples  Number  of  samples  in  input  vectors 

bins  Number  of  bins  in  the  Probability  Density  Function  (PDF) 

samples _per_symbol  Number  of  samples  per  data  symbol 

numjrefs  Number  of  reference  signals 

m_type  Modulation  type 
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Outputs: 

mnt 

g 

out 

clk_out 

hold_vec 

pdf 

bp 


Plot  of  the  Memoryless  Nonlinear  Transform 
Complex  transformed  sequence 
Decision  of  nonlinearity 

Timing  signal:  goes  high  when  valid  data  sample  is  available  at  the  output 
of  cpa/s_cdf3_nonlin 

Timing  signal:  goes  low  when  vector  outputs  are  available 

Plot  of  the  PDF  of  magnitude 

Breakpoints  (PDF  and  MNT  interval  boundaries) 


See  also: 

cpa/cdO_nonlin,  histo/correlator2,  romelib/timing,  rl/inf_vsource 


Name: 

cpa/s_pdf_gtjil_nct 

Description: 

This  hierarchical  block  buffers  the  received  data  and  converts  it  from  serial  to  vector  data  and 
applies  it  to  the  cpa/pdf_gt_nonlin  nonlinearity.  This  block  does  not  generate  the  summation 
part  of  the  LO  decision  statistic. 

The  romelib/timing  block  generates  timing  waveforms  which  enable  and  disable  the 
nonlinearity  and  all  subsequent  blocks  at  appropriate  times.  Of  particular  interest  are  the 
CLK_OUT  and  HOLD_VEC  signals.  The  HOLD_VEC  signal  holds  vector  blocks  while  the 
serial  data  is  buffered  into  vector  form.  This  signal  goes  low  during  the  simulation  iteration  in 
which  the  vector  processing  completes  and  the  outputs  are  available.  The  CLK_OUT  signal 
provides  flexibility  for  use  in  multirate  systems  by  pulsing  high  when  a  valid  data  sample  is 
produced  at  the  output  of  the  cpa/sjxlf _gt_nl_nc  block.  The  CLK_OUT  signal  has  the  same 
period  as  the  CLK_IN  signal,  thus  output  samples  are  produced  at  the  same  rate  as  that  of  the 
input  signal.  The  remaining  romelib/timing  signals  are  internal  to  the  cpa/s_pdf_gt_nl_nc  block, 
and  are  used  in  serial  to  vector  and  vector  to  serial  buffering. 

The  SAMPLES_PER_SYMBOL,  NUM_REFS,  and  M_TYPE  are  characteristics  of  the 
reference  symbols.  The  reference  signals  must  be  created  before  using  this  block,  and  are  stored 
in  the  /spwdata/ref  library.  This  permits  the  designer  to  use  any  desired  signal  constellation. 
The  NUM_REFS  parameter  must  be  exported  to  the  histo/correlator2  block,  and  the  other  two 
parameters  are  exported  to  the  filename  of  the  rl/inf_vsource  block  as  "ref/"&substr 
C  m_typeH :  model  ,2,length("  m_type" :  model)-2)&"_H  &xstring("  samples_per_symbol" :  model)& 
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"sPb_i’"  for  the  I  channel  signals  and  “  ref/ " &substr( "  m_type “ : model , 2 , length( * mtype " : 
model)-2)  &’,_"&xstring(Msamples_per_symbor:model)&"spb_q’"  for  the  Q  channel  signals. 
For  example,  the  Filenames  for  a  Binary  Phase  Shift  Keying  (BPSK)  signal  with  25  samples  per 
symbol  would  be  /spwdata/ref/bpsk_25spb_i  and  /spwdata/ref/bpsk_25spb_q. 

Inputs: 

in  Input  message  sequence  (complex  input) 

clk_in  Clock  input  for  romelib/timing 

Parameters: 

samples  Number  of  samples  in  input  vectors 

bins  Number  of  bins  in  the  Probability  Density  Function  (PDF) 

samples_per_symboi  Number  of  samples  per  data  symbol 

num_refs  Number  of  reference  signals 

mjype  Modulation  type 

Outputs: 
mnt 
g 

out 

clk_out 

hold_vec 
pdf 
bp 


See  also: 

cpa/pdf_gt_nonlin,  histo/correlator2,  romelib/timing,  cpa/s_pdf_gt_nonli,  rl/inf_vsource 


Plot  of  the  Memoryless  Nonlinear  Transform 
Complex  transformed  sequence  (vector) 

Complex  transformed  sequence  (serial) 

Timing  signal:  goes  high  when  valid  data  sample  is  available  at  the  output 
of  cpa/s_pdf_gt_nl_nc 

Timing  signal:  goes  low  when  vector  outputs  are  available 

Plot  of  the  PDF  of  magnitude 

Breakpoints  (PDF  and  MNT  interval  boundaries) 


Name: 

cpa/  s_pd  f__gt_n°nlit 
Description: 

This  hierarchical  block  buffers  the  received  data  and  converts  it  from  serial  to  vector  data  and 
applies  it  to  the  cpa/pdf_gt_nonlin  nonlinearity.  In  addition,  the  summation  in  the  LO  decision 
statistic  is  formed  using  the  histo/correlator2  block.  The  output  of  the  correlator  is  converted 
back  to  serial  form. 


The  romelib/ timing  block  generates  timing  waveforms  which  enable  and  disable  the 
nonlinearity  and  all  subsequent  blocks  at  appropriate  times.  Of  particular  interest  are  the 
CLK_OUT  and  HOLD_VEC  signals.  The  HOLD_VEC  signal  holds  vector  blocks  while  the 
serial  data  is  buffered  into  vector  form.  This  signal  goes  low  during  the  simulation  iteration  in 
which  the  vector  processing  completes  and  the  outputs  are  available.  The  CLK_OUT  signal 
provides  flexibility  for  use  in  multirate  systems  by  pulsing  high  when  a  valid  data  sample  is 
produced  at  the  output  of  the  cpa/sjpdf_gt_nonli  block.  The  CLK_OUT  signal  has  the  same 
period  as  the  CLK_IN  signal,  thus  output  samples  are  produced  at  the  same  rate  as  that  of  the 
input  signal.  The  remaining  romelib/timing  signals  are  internal  to  the  cpa/s_pdf_gt_nonli  block, 
and  are  used  in  serial  to  vector  and  vector  to  serial  buffering. 

The  SAMPLES_PER_SYMBOL,  NUM_REFS,  and  M  TYPE  are  characteristics  of  the 
reference  symbols.  The  reference  signals  must  be  created  before  using  this  block,  and  are  stored 
in  the  /spwdata/ref  library.  This  permits  the  designer  to  use  any  desired  signal  constellation. 
The  NUM_REFS  parameter  must  be  exported  to  the  histo/correlator2  block,  and  the  other  two 
parameters  are  exported  to  the  filename  of  the  rl/inf_vsource  block  as  "ref/"&substr 
( "  m_type " :  model  ,2 ,  length  ("  m_type  ” :  model)-2)&"_''  &xstring( "  samples_per_symbol" :  model)& 
Mspb_i’"  for  the  I  channel  signals  and  " ref/ " &substr( * m_type " : model , 2 , length( M m_type" : 
model)-2)  &"_"&xstring("samples_per_symbor:model)&"spb_q,M  for  the  Q  channel  signals. 
For  example,  the  filenames  for  a  Binary  Phase  Shift  Keying  (BPSK)  signal  with  25  samples  per 
symbol  would  be  /spwdata/ref/bpsk_25spb_i  and  /spwdata/ref/bpsk_25spb_q. 


Inputs: 

in  Input  message  sequence  (complex  input) 

clk_in  Clock  input  for  romelib/timing 


Parameters: 

samples 

bins 

samples_per_symbol 

numjrefs 

m_type 


Number  of  samples  in  input  vectors 

Number  of  bins  in  the  Probability  Density  Function  (PDF) 

Number  of  samples  per  data  symbol 

Number  of  reference  signals 

Modulation  type 
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Outputs: 

mnt 

g 

out 

elkjMIt 

hold_vec 

pdf 

bp 


Plot  of  the  Memoryless  Nonlinear  Transform 
Complex  transformed  sequence 
Decision  of  nonlinearity 

Timing  signal:  goes  high  when  valid  data  sample  is  available  at  the  output 
of  cpa/sjxlf _gt_nonli 

Timing  signal:  goes  low  when  vector  outputs  are  available 

Plot  of  the  PDF  of  magnitude 

Breakpoints  (PDF  and  MNT  interval  boundaries) 


See  also: 

cpa/pdf _gt_nonlin,  histo/correlator2,  romelib/timing,  rl/inf_v$ource 


Name: 

cpa/slopel 

Description: 

This  block  computes  the  three  point  derivative  of  the  histogram  CDF  for  use  in  the  CPA 
MNT  construction.  Extra  bins  are  added  to  compute  the  slopes  at  the  endpoints.  This  is  the 
method  described  in  Appendix  3  of  IIT  Final  Report  1991,  F30602-91-C-0059.  This  block  is 
used  in  the  cpa/cdf2  and  cpa/cdf3  blocks. 

Note: 

The  additional  derivative  at  each  end  is  always  zero,  since 
HISTO_CDF[-2] = HISTO_CDF[OJ =0.0  and  HISTO_CDF[BINS]  =HISTO_CDF[BINS  +2] = 1 . 0. 

Inputs: 

bp  Breakpoints  (CDF  interval  boundaries) 

histo_cdf  Histogram  CDF 

Parameters: 

bins  Number  of  bins 

Outputs: 
df 
ddf 


First  derivative 
Second  derivative 


Name: 

cpa/slope2 

Description: 

This  block  computes  the  three  point  derivative  of  the  histogram  PDF  for  use  in  the  CPA  MNT 
construction.  Extra  bins  are  added  to  compute  the  slopes  at  the  endpoints.  This  block  is  used 
in  the  cpa/pdf_gt  block. 

Note: 

The  additional  derivative  at  each  end  is  always  zero,  since 
HISTO_PDF[-2] = HISTO_PDF[0]  =O.QmdHISTO_PDF[BINS] = HISTO_PDF[BINS  +2]  =0.0. 

Inputs: 

bp  Breakpoints  (PDF  interval  boundaries) 

histojxlf  Histogram  PDF 

Parameters: 

bins  Number  of  bins 

Outputs: 
df 
ddf 


First  derivative 
Second  derivative 


Name: 

fourier/mntf 

Description: 

This  block  obtains  the  Memoryless  Nonlinear  Transform  (MNT)  from  a  Fourier  Series 
Approximation  (FSA)  of  a  Probability  Density  Function  (PDF).  It  is  a  hierarchical  block  which 
contains  two  fourier/mnt_calc  blocks.  One  of  them  computes  the  MNT  for  each  point  of  the 
magnitude  of  the  input  message  sequence,  and  the  other  generates  a  plot  of  the  MNT  from  the 
minimum  to  the  maximum  magnitude  value  with  the  aid  of  a  vec/minmax_ramp  block. 

This  block  operates  on  the  magnitude  of  the  two  dimensional  message  signal. 

Inputs: 

a  Vector  of  A  coefficients  (cosine  terms) 

b  Vector  of  B  coefficients  (sine  terms) 

T  "Period-  of  the  FSA,  equal  to  MAX  -  MIN  (T  is  computed  by  the 

fourier/pdf  block) 

message  Input  message  sequence  (2  dimensional) 

Parameters: 

samples  Number  of  samples  in  data  vector 

points  Number  of  plot  points  for  MNT 

P  Order  of  FSA 

Outputs: 

g_val  MNT  of  each  data  sample 

mnt  Plot  of  MNT  from  minimum  to  maximum  data  value 

See  also: 

fourier/mnt_calc,  fourier/pdf,  vec/minmax_ramp 


Name: 

fourier/mnt_calc 

Description: 

This  block  obtains  the  Memoryless  Nonlinear  Transform  (MNT)  from  a  Fourier  Series 
Approximation  (FSA)  of  a  Probability  Density  Function  (PDF).  The  MNT  of  each  point  of  the 
input  message  sequence  is  stored  in  the  G_VAL  output. 
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This  block  operates  on  the  magnitude  of  the  two  dimensional  message  signal. 

Inputs: 
a 
b 
T 

message 
Parameters: 

samples  Number  of  samples  in  data  vector 

P  Order  of  FSA 

Outputs: 

g_val  MNT  of  each  data  sample 


Name: 

fourier/nonlin+ 

Description: 

This  is  a  hierarchical  block  which  applies  the  nonlinearity  to  the  input  signal  in  vector  format. 
The  histo/rec_to_polar  block  converts  the  two  dimensional  input  signal  from  rectangular  to  polar 
coordinates,  and  a  histogram  approximation  of  the  Probability  Density  Function  (PDF)  of  the 
magnitude  is  computed  by  the  histo/histo  block.  The  fourier/pdf  block  then  constructs  a  Fourier 
Series  Approximation  (FSA)  of  the  PDF  based  on  the  histogram.  The  fourier/mnt  block  obtains 
the  Memoryless  Nonlinear  Transform  (MNT)  from  the  FSA  PDF  and  applies  it  to  the  magnitude 
of  the  input  message  sequence.  The  result  is  converted  back  to  rectangular  coordinates  by  the 
histo/polar_to_rec  block  using  the  unmodified  phase.  The  vec/ramp  block  computes  the 
breakpoints  for  plotting  purposes. 

Refer  to  the  fourier/mnt  block  for  more  information  on  the  MNT. 

Inputs: 

i  Jn  In-Phase  component  of  input  message  sequence 

qjn  Quadrature  component  of  input  message  sequence 

Parameters: 

samples  Number  of  samples  in  input  vectors 

bins  Number  of  bins  in  the  histogram  PDF 


Vector  of  A  coefficients  (cosine  terms) 

Vector  of  B  coefficients  (sine  terms) 

"Period"  of  the  FSA,  equal  to  MAX  -  MIN  (T  is  computed  by  the 
fourier/pdf  block) 

Input  message  sequence  (2  dimensional) 


Outputs: 

pdf  Plot  of  the  FSA  PDF 

i_out  In-Phase  component  of  transformed  sequence 

q_out  Quadrature  component  of  transformed  sequence 

mnt  Plot  of  the  FSA  MNT 

bp  Breakpoints  (PDF  and  MNT  interval  boundaries) 

See  also: 

fourier/pdf,  fourier/pdfjplot,  fourier/mnt,  histo/histo,  histo/polar_to_rec, 
histo/rec_to_polar,  vec/ramp 


Name: 

fourier/pdf 

Description: 

This  block  obtains  the  Fourier  Series  Approximation  (FSA)  of  a  Probability  Density  Function 
(PDF).  This  FSA  method  is  based  on  a  histogram  of  the  input  data  magnitude. 


Inputs: 

bins  Histogram  PDF 

min  Minimum  data  value 

max  Maximum  data  value 


Parameters: 

samples  Number  of  samples  in  data  vector 

P  Order  of  FSA 


Outputs: 

a  Vector  of  A  coefficients  (cosine  terms) 

b  Vector  of  B  coefficients  (sine  terms) 

T  "Period"  of  the  FSA,  equal  to  MAX  -  MIN 
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Name: 

fourier/pdf_plott 

Description: 

This  hierarchical  block  plots  the  Fourier  Series  Approximation  (FSA)  of  a  Probability  Density 
Function  (PDF).  A  vec/minmax_ramp  block  generates  a  ramp  from  the  minimum  to  the 
maximum  input  values,  and  the  fourier/pdfplot  block  computes  the  PDF  approximation  for  each 
point  in  the  ramp. 

Inputs: 
a 
b 
T 

message 

Parameters: 
samples 
points 
P 

Outputs: 
pdf 

See  also: 

fourier/pdfplot,  vec/minmax_ramp 


Vector  of  A  coefficients  (cosine  terms) 

Vector  of  B  coefficients  (sine  terms) 

"Period"  of  the  FSA,  equal  to  MAX  -  MIN  (T  is  computed  by  the 
fourier/pdf  block) 

Input  message  sequence 


Number  of  samples  in  data  vector 
Number  of  plot  points  for  PDF 
Order  of  FSA 


FSA  PDF 


Name: 

fourier/pdfplot 

Description: 

This  block  plots  the  Fourier  Series  Approximation  (FSA)  of  a  Probability  Density  Function 
(PDF).  The  PDF  approximation  is  computed  for  each  point  in  the  input,  which  must  be  a  linear 
ramp  over  the  desired  plot  range. 
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Inputs: 

points  Input  data  vector 

a  Vector  of  A  coefficients  (cosine  terms) 

b  Vector  of  B  coefficients  (sine  terms) 

T  "Period"  of  the  FSA,  equal  to  MAX  -  MIN  (T  is  computed  by  the 

fourier/pdf  block) 


Parameters: 

points  Number  of  plot  points  for  PDF 

P  Order  of  FSA 


Outputs: 

pdf  FSA  PDF 


Name: 

fourier/s_nonlin_nct 

Description: 

This  hierarchical  block  buffers  the  received  data  and  converts  it  from  serial  to  vector  data  and 
applies  it  to  the  fourier/nonlin  nonlinearity.  This  block  does  not  generate  the  summation  part 
of  the  LO  decision  statistic. 

The  romelib/timing  block  generates  timing  waveforms  which  enable  and  disable  the 
nonlinearity  and  all  subsequent  blocks  at  appropriate  times.  Of  particular  interest  are  the 
CLK_OUT  and  HOLD_VEC  signals.  The  HOLD_VEC  signal  holds  vector  blocks  while  the 
serial  data  is  buffered  into  vector  form.  This  signal  goes  low  during  the  simulation  iteration  in 
which  the  vector  processing  completes  and  the  outputs  are  available.  The  CLK_OUT  signal 
provides  flexibility  for  use  in  multirate  systems  by  pulsing  high  when  a  valid  data  sample  is 
produced  at  the  output  of  the  fourier/s_nonlin_nc  block.  The  CLK_OUT  signal  has  the  same 
period  as  the  CLK_IN  signal,  thus  output  samples  are  produced  at  the  same  rate  as  that  of  the 
input  signal.  The  remaining  romelib/timing  signals  are  internal  to  the  fourier/s_nonlin_nc  block, 
and  are  used  in  serial  to  vector  and  vector  to  serial  buffering. 

The  SAMPLES_PER_SYMBOL,  NUM_REFS,  and  M_TYPE  are  characteristics  of  the 
reference  symbols.  The  reference  signals  must  be  created  before  using  this  block,  and  are  stored 
in  the  /spwdata/ref  library.  This  permits  the  designer  to  use  any  desired  signal  constellation. 
The  NUM_REFS  parameter  must  be  exported  to  the  histo/corre!ator2  block,  and  the  other  two 
parameters  are  exported  to  the  filename  of  the  rl/inf_vsource  block  as  "ref/"&substr 
(M  mjype" :  model  ,2,length("m_type” :  model)-2)&"_"&xstring("samples_per_symbor :  model)& 


A-24 


HspbJ’-  for  the  I  channel  signals  and  " ref/ "  &  substr  ( “ m_type " : model ,2 , length( " m_type " : 
model)-2)  &"_"&xstring(,,samples_per_symbor:model)&"spb_q’"  for  the  Q  channel  signals. 
For  example,  the  filenames  for  a  Binary  Phase  Shift  Keying  (BPSK)  signal  with  25  samples  per 
symbol  would  be  /spwdata/ref/bpsk_25spb_i  and  /spwdata/ref/bpsk_25spb_q. 

Inputs: 

in  Input  message  sequence  (complex  input) 

clk_in  Clock  input  for  romelib/timing 

Parameters: 
samples 
bins 

samples_per_symbol 
num_refs 
m_type 

Outputs: 

mnt  Plot  of  the  Memoryless  Nonlinear  Transform 

g  Complex  transformed  sequence  (vector) 

out  Complex  transformed  sequence  (serial) 

clk_out  Timing  signal:  goes  high  when  valid  data  sample  is  available  at  the  output 

of  fourier/s_nonlin_nc 

hold_vec  Timing  signal:  goes  low  when  vector  outputs  are  available 

pdf  Plot  of  the  PDF  of  magnitude 

bp  Breakpoints  (PDF  and  MNT  interval  boundaries) 


Number  of  samples  in  input  vectors 
Number  of  bins  in  the  histogram 
Number  of  samples  per  data  symbol 
Number  of  reference  signals 
Modulation  type 


See  also: 

fourier/nonlin,  romelib/timing,  rl/inf_vsource,  fourier/serial_nonlin 
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Name: 

fourier/serial_nonlint 

Description: 

This  hierarchical  block  buffers  the  received  data  and  converts  it  from  serial  to  vector  data  and 
applies  it  to  the  fourier/nonlin  nonlinearity.  In  addition,  the  summation  in  the  LO  decision 
statistic  is  formed  using  the  histo/correlator2  block.  The  output  of  the  correlator  is  converted 
back  to  serial  form. 


The  romelib/timing  block  generates  timing  waveforms  which  enable  and  disable  the 
nonlinearity  and  all  subsequent  blocks  at  appropriate  times.  Of  particular  interest  are  the 
CLK_OUT  and  HOLD_VEC  signals.  The  HOLD_VEC  signal  holds  vector  blocks  while  the 
serial  data  is  buffered  into  vector  form.  This  signal  goes  low  during  the  simulation  iteration  in 
which  the  vector  processing  completes  and  the  outputs  are  available.  The  CLK_OUT  signal 
provides  flexibility  for  use  in  multirate  systems  by  pulsing  high  when  a  valid  data  sample  is 
produced  at  the  output  of  the  fourier/serial  jionlin  block.  The  CLK_OUT  signal  has  the  same 
period  as  the  CLK_IN  signal,  thus  output  samples  are  produced  at  the  same  rate  as  that  of  the 
input  signal.  The  remaining  romelib/timing  signals  are  internal  to  the  fourier/serial  jionlin 
block,  and  are  used  in  serial  to  vector  and  vector  to  serial  buffering. 

The  SAMPLES_PER_SYMBOL,  NUM_REFS,  and  M_TYPE  are  characteristics  of  the 
reference  symbols.  The  reference  signals  must  be  created  before  using  this  block,  and  are  stored 
in  the  /spwdata/ref  library.  This  permits  the  designer  to  use  any  desired  signal  constellation. 
The  NUM_REFS  parameter  must  be  exported  to  the  histo/correlator2  block,  and  the  other  two 
parameters  are  exported  to  the  filename  of  the  rl/inf_vsource  block  as  "ref/"&substr 
( "  m_type " :  model ,  2 ,  length("  m_type" :  model)-2)&  "_"&xstring("  samples _per_symbol " :  model)& 
”spb_i’"  for  the  I  channel  signals  and  ,,ref/"&substr("m_type,,:model,2,length("m_type": 
model)-2)  &"_"&xstring("samplesjper_symbor:modei)&"spb_q’"  for  the  Q  channel  signals. 
For  example,  the  filenames  for  a  Binary  Phase  Shift  Keying  (BPSK)  signal  with  25  samples  per 
symbol  would  be  /spwdata/ref/bpsk_25spb_i  and  /spwdata/ref/bpsk_25spb_q. 


Inputs: 

in  Input  message  sequence  (complex  input) 

clk_in  Clock  input  for  romelib/timing 


Parameters: 

samples 

bins 

samples_per_symbol 

num_refs 

m_type 


Number  of  samples  in  input  vectors 
Number  of  bins  in  the  histogram 
Number  of  samples  per  data  symbol 
Number  of  reference  signals 
Modulation  type 
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Outputs: 

mnt 

g 

out 

clk_out 

hold_vec 

pdf 

bp 


See  also: 

fourier/nonlin. 


Plot  of  the  Memoryless  Nonlinear  Transform 
Complex  transformed  sequence 
Decision  of  nonlinearity 

Timing  signal:  goes  high  when  valid  data  sample  is  available  at  the  output 
of  fourier/serialjionlin 

Timing  signal:  goes  low  when  vector  outputs  are  available 

Plot  of  the  PDF  of  magnitude 

Breakpoints  (PDF  and  MNT  interval  boundaries) 


histo/correlator2,  romelib/ timing,  rl/inf_vsource 
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Name: 

histo/ave_bp 

Description: 

The  histogram  is  defined  only  between  breakpoints,  leaving  values  at  the  breakpoints 
undefined.  However,  some  of  the  polynomial  curve  fitting  algorithms  require  knowledge  of  the 
histogram  values  at  the  breakpoints.  This  function  computes  the  values  at  the  breakpoints  by 
averaging  the  values  of  the  adjacent  bins.  The  first  and  last  values  are  taken  to  be  the  values 
of  the  first  and  last  bins.  The  output  "histogram"  will  have  one  more  data  point  than  the  input 
histogram. 

Inputs: 

in  Input  histogram 

Parameters: 

bins  Number  of  bins 

Outputs: 

out  Output  "histogram" 


Name: 

histo/correnf  )r2 
Description: 

This  block  implements  the  matched  filter  -  correlator  for  a  two  dimensional  system.  The  input 
reference  signals,  I_REF  and  Q_REF,  are  assumed  to  have  all  possible  signal  pairs  stored 
sequentially.  The  width  of  each  signal  is  given  by  SAMPLES_PER_SYMBOL.  The  total 
number  of  reference  signals  is  given  by  NUM_REFS. 

Inputs: 
i 

q 

i_ref 
q_ref 

Parameters: 

samples  Number  of  samples  in  input  vectors 

sampies_per_symbol  Number  of  samples  per  data  symbol 
num__refs  Number  of  reference  symbols 


In-Phase  component  of  input  message  sequence 
Quadrature  component  of  input  message  sequence 
I-channel  Reference  symbols 
Q-channel  Reference  symbols 
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Outputs: 

i_decision  Decision  of  In-Phase  component  of  the  message 

q_decision  Decision  of  Quadrature  component  of  the  message 


Name: 

histo/enonlin* 

Description: 

This  is  a  hierarchical  block  which  applies  the  nonlinearity  to  the  input  signal  in  vector  format. 
The  histo/rec_to _polar  block  converts  the  two  dimensional  input  signal  from  rectangular  to  polar 
coordinates,  and  the  data  is  sorted  by  the  vec/heap_sort  block.  An  equiprobable  bin  histogram 
approximation  of  the  Probability  Density  Function  (PDF)  of  the  magnitude  is  computed  by  the 
histo/equi  block.  The  histo/mnt2  block  obtains  the  Memoryless  Nonlinear  Transform  (MNT) 
from  the  histogram  and  applies  it  to  the  magnitude  of  the  input  message  sequence.  The  result 
is  convened  back  to  rectangular  coordinates  by  the  histo/polar_to_rec  block  using  the  unmodified 
phase. 

Refer  to  the  histo/mnt2  block  for  more  information  on  the  MNT. 

Inputs: 

i_in  In-Phase  component  of  input  message  sequence 

q_in  Quadrature  component  of  input  message  sequence 

Parameters: 

samples  Number  of  samples  in  input  vectors 

bins  Number  of  bins  in  the  histogram  PDF 

Outputs: 
pdf 
i_out 
q_out 
mnt 
bp 

See  also: 

vec/heap_sort,  histo/equi,  histo/mnt2,  hi$to/rec_to__polar,  histo/polar_to_rec 


Equiprobable  bin  histogram  PDF 
In-Phase  component  of  transformed  sequence 
Quadrature  component  of  transformed  sequence 
Plot  of  the  Memoryless  Nonlinear  Transform 
Breakpoints  (PDF  and  MNT  interval  boundaries) 
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Name: 

histo/equi 

Description: 

This  block  generates  an  equiprobable  histogram  Probability  Density  Function  (PDF)  and  an 
equiprobable  histogram  Cumulative  Distribution  Function  (CDF)  of  the  input  data  sequence. 
Since  The  PDF  and  CDF  are  equiprobable,  the  probability  of  a  data  point  falling  in  any  bin  is 
equal  to  1/BINS.  THE  INPUT  DATA  MUST  BE  SORTED  IN  ASCENDING  ORDER 
BEFORE  IT  IS  APPLIED  TO  THIS  BLOCK. 

Inputs: 
data 

Parameters: 
bins 
samples 

Outputs: 
bp 
pdf 
cdf 


Input  data  sequence  (must  be  sorted  in  ascending  order) 


Number  of  bins  to  generate 
Number  of  samples  in  data  vector 


Breakpoints  (PDF  and  CDF  interval  boundaries) 
Equi-probable  bin  histogram  PDF 
Equi-probable  bin  histogram  CDF 


Name: 

histo/histo 

Description: 

This  block  generates  a  histogram  Probability  Density  Function  (PDF)  of  the  input  data 
sequence.  The  width  of  each  bin  is  (MAX-MIN)/BINS.  The  output  is  scaled  so  that  the  total 
area  of  the  histogram  is  equal  to  1. 

Inputs: 

message  Input  data  sequence 

Parameters: 

message  Jen  Number  of  samples  in  data  vector 

bins  Number  of  bins  to  generate 
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Outputs: 

bins  Probability  Density  Function 

min  Minimum  data  value 

max  Maximum  data  value 


Name: 

histo/mnt2f 

Description: 

This  block  obtains  the  Memoryless  Nonlinear  Transform  (MNT)  from  a  histogram  Probability 
Density  Function  (PDF).  It  is  a  hierarchical  block  which  contains  a  histo/mnt_calc2  block  and 
a  histo/mnt_out2  block.  The  mnt_calc2  block  computes  the  MNT,  and  the  mnt_out2  block 
applies  the  MNT  to  each  data  point  of  the  magnitude  of  the  input  message  sequence. 

This  block  is  a  generalized  version  of  histo/mntjq;  it  has  a  breakpoint  input  so  it  may  be  used 
for  a  histogram  with  bins  of  arbitrary  width.  Thus,  it  may  be  used  in  conjunction  with  either 
a  histogram  or  an  equiprobable  bin  histogram.  The  MNT  of  the  K*  histogram  bin  is 

rfn ,  1  .  log(PDF[K+l])-log(PDF[K-l]) 

J  BP[K]  BP[K]-BP[K-1] 

In  the  case  of  a  one  dimensional  data  signal,  the  1/BP[K]  term  is  not  included  in  g[K].  For 
a  two  dimensional  data  signal,  this  block  operates  only  on  the  magnitude. 

The  G_PLOT  output  is  the  MNT  computed  for  each  histogram  bin.  Thus,  the  length  of  the 
G_PLOT  vector  must  be  the  same  as  the  BINS  vector.  The  length  of  the  G_VAL  vector  is  the 
same  as  the  DATA  vector. 


Inputs: 

bins  Histogram  PDF 

bp  Breakpoints  (PDF  and  MNT  interval  boundaries) 

data  Input  message  sequence  (1  or  2  dimensional) 


Parameters: 

bins  Number  of  bins  in  the  PDF 

samples  Number  of  samples  in  data  vector 

dim  Dimensionality  of  the  message  sequence  (1  or  2) 


A-31 


Outputs: 

g_val  MNT  of  each  data  sample 

g_plot  Plot  of  MNT  from  minimum  to  maximum  data  value 

See  also: 

histo/mnt_calc2,  histo/mnt_out2 


Name: 

histo/mnt_calc2 

Description: 

This  function  computes  the  Memoryless  Nonlinear  Transform  (MNT)for  a  histogram  with  bins 
of  arbitrary  width.  The  MNT  of  the  K*  histogram  bin  is 

cfK1  -  1  ,  log(PDF[K+l])-log(PDF[K-l]) 

61  BP[K]  BP[K]-BP[K-1] 

In  the  case  of  a  one  dimensional  data  signal,  the  1/BP[KJ  term  is  not  included  in  g[K].  For 
a  two  dimensional  data  signal,  this  block  operates  only  on  the  magnitude. 


Inputs: 

bins  Histogram  Probability  Density  Function  (PDF) 

bp  Breakpoints  (PDF  and  MNT  interval  boundaries) 


Parameters: 

bins  Number  of  bins  in  the  PDF 

dim  Dimensionality  of  the  message  sequence  (1  or  2) 


Outputs: 

g 


MNT  for  each  histogram  bin 
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Name: 

histo/mnt_calc_2d 

Description: 

This  function  computes  the  Memoryless  Nonlinear  Transform  (MNT)  for  a  histogram  with 
bins  of  equal  width.  The  MNT  of  the  K*  histogram  bin  is 

gfKl3J_  log(PDF[K+l])-log(PDF[K-l]) 

6  R  WIDTH 

where  WIDTH  =  MAX  -  MIN  and 
R  *  MIN  +  K  *  WIDTH. 

This  block  operates  on  the  magnitude  of  the  two  dimensional  data  signal. 


Inputs: 

bins  Histogram  Probability  Density  Function  (PDF) 

min  Minimum  magnitude  of  input  message  sequence 

max  Maximum  magnitude  of  input  message  sequence 


Parameters: 

bins  Number  of  bins  in  the  PDF 


Outputs: 

g 


MNT  for  each  histogram  bin 


Name: 

histo/mnt_iqt 

Description: 

This  block  obtains  the  Memoryless  Nonlinear  Transform  (MNT)  from  a  histogram  Probability 
Density  Function  (PDF).  It  is  a  hierarchical  block  which  contains  a  histo/mnt_calc_d2  block 
and  a  histo/mnt_out  block.  The  mnt_calc_2d  block  computes  the  MNT,  and  the  mnt_out  block 
applies  the  MNT  to  each  data  point  of  the  magnitude  of  the  input  message  sequence. 

This  block  assumes  that  each  histogram  bin  has  equal  width.  In  addition,  input  message 
sequence  is  assumed  to  be  the  magnitude  of  a  two  dimensional  signal.  Tl*#  MNT  of  the  K* 
histogram  bin  is 
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J_  log(PDF[K+l])-log(PDF[K-l]) 

151  J  R  WIDTH 

where  WIDTH  =  MAX  -  MIN  and 
R  *  MIN  +  K  *  WIDTH. 

The  G_PLOT  output  is  the  MNT  computed  for  each  histogram  bin.  Thus,  the  length  of  the 
G_PLOT  vector  must  be  the  same  as  the  BINS  vector.  The  length  of  the  G_VAL  vector  is  the 
same  as  the  MESSAGE  vector. 


Inputs: 

bins  Histogram  PDF 

min  Minimum  magnitude  of  input  message  sequence 

max  Maximum  magnitude  of  input  message  sequence 

message  Input  message  sequence  (2  dimensional) 


Parameters: 

bins  Number  of  bins  in  the  PDF 

samples  Number  of  samples  in  data  vector 


Outputs: 

g_val  MNT  of  each  data  sample 

g_plot  Plot  of  MNT  from  minimum  to  maximum  data  value 

See  also: 

histo/mnt_calc_2d,  histo/mnt_out 


Name: 

histo/mnt_out 

Description: 

This  block  applies  the  Memoryless  Nonlinear  Transform  (MNT)  generated  by  the 
histo/mnt_calc_2d  block  to  each  sample  in  the  input  message  sequence.  This  block  assumes  that 
each  histogram  bin  has  equal  width.  In  addition,  the  input  message  sequence  is  assumed  to  be 
the  magnitude  of  a  two  dimensional  signal. 
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Inputs: 

g  MNT  for  input  message  sequence 

min  Minimum  magnitude  of  input  message  sequence 

max  Maximum  magnitude  of  input  message  sequence 

message  Input  message  sequence  (2  dimensional) 

Parameters: 

bins  Number  of  bins  in  the  Probability  Density  Function 

samples  Number  of  samples  in  data  vector 

Outputs: 

g_val  MNT  of  each  input  sample 


Name: 

histo/mnt_out2 

Description: 

This  block  applies  the  Memoryless  Nonlinear  Transform  (MNT)  generated  by  the 
histo/mnt_calc2  block  to  each  sample  in  the  input  message  sequence.This  block  permits  the  use 
of  arbitrarily  spaced  breakpoints. 

Inputs: 

a 
© 

bp 
data 

Parameters: 

bins  Number  of  bins  in  the  histogram 

samples  Number  of  samples  in  data  vector 

Outputs: 

g_val  MNT  of  each  data  sample 


MNT  for  input  message  sequence 
Breakpoints  (MNT  interval  boundaries) 

Input  message  sequence  (1  or  2  dimensional) 
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Name: 

histo/nonlin6f 

Description: 

This  is  a  hierarchical  block  which  applies  the  nonlinearity  to  the  input  signal  in  vector  format. 
The  histo/rec_to_polar  block  converts  the  two  dimensional  input  signal  from  rectangular  to  polar 
coordinates,  and  a  histogram  approximation  of  the  Probability  Density  Function  (PDF)  of  the 
magnitude  is  computed  by  the  histo/histo  block.  The  histo/mntjq  block  obtains  the  Memoryless 
Nonlinear  Transform  (MNT)  from  the  histogram  and  applies  it  to  the  magnitude  of  the  input 
message  sequence.  The  result  is  converted  back  to  rectangular  coordinates  by  the 
histo/polar_to_rec  block  using  the  unmodified  phase.  The  vec/ramp  block  computes  the 
breakpoints  for  plotting  purposes. 

Refer  to  the  histo/mntjq  block  for  more  information  on  the  MNT. 

Inputs: 
ijn 
q_in 

Parameters: 
samples 
bins 

Outputs: 
pdf 
i_out 
q_out 
mnt 
bp 

See  also: 

histo/histo,  histo/mntjq,  histo/recjo_polar,  histo/polar_to_rec,  vec/ramp 


In-Phase  component  of  input  message  sequence 
Quadrature  component  of  input  message  sequence 


Number  of  samples  in  input  vectors 
Number  of  bins  in  the  histogram  PDF 


Histogram  PDF 

In-Phase  component  of  transformed  sequence 
Quadrature  component  of  transformed  sequence 
Plot  of  the  Memoryless  Nonlinear  Transform 
Breakpoints  (PDF  and  MNT  interval  boundaries) 
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Name: 

histo/polar_to_rec 

Description: 

This  block  accepts  a  complex  vector  in  polar  coordinates  and  converts  it  to  a  complex  vector 
in  rectangular  coordinates  as  follows: 

I=R*cos(A) 

Q=R*sin(A) 

where  R  is  the  magnitude  and  A  is  the  phase. 


Inputs: 

mag  Magnitude  of  input  data 

phase  Phase  of  input  data 


Parameters: 

samples  Number  of  samples  in  input  and  output  vectors 

Outputs: 

i  In-Phase  component  of  input  data 

q  Quadrature  component  of  input  data 


Name: 

histo/prob_errort 

Description: 

This  hierarchical  block  computes  the  Probability  of  Bit  Error  (Pb)  of  a  serial  data  stream. 
The  number  of  errors  in  both  the  I  and  Q  channel  are  summed  over  SYMBOLS  JPER_CALC 
data  bits  in  each  channel.  The  Pb  is  equal  to  this  sum  divided  by  (SYMBOLS_PER_CALC*2). 

The  CLOCKJN  control  signal  must  go  high  when  a  sample  is  present  at  the  input.  This 
signal  may  be  generated  using  the  CLK_OUT  output  of  the  romelib/timing  block. 
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Inputs: 

ijn 

i_dec 

q_dec 

q_in 

clockin 


Reference  signal— I  channel 
Decision  of  receiver--I  channel 
Decision  of  receiver--Q  channel 
Reference  signal~Q  channel 
Clock  enable 


Parameters: 

symbols_per_calc  Number  of  symbols  per  Pb  calculation 
sampies_per_symbol  Number  of  samples  per  data  symbol 


Outputs: 

pe  Probability  of  Bit  Error 

pe_clock  Clock  signal  for  Pb  signal  sink 


Name: 

histo/rec_to_polar 

Description: 

This  block  accepts  a  complex  vector  in  rectangular  coordinates  and  converts  it  to  a  complex 
vector  in  polar  coordinates  as  follows: 

R=sqrt(IA2+QA2) 

A=atan(Q/I) 

where  R  is  the  magnitude  and  A  is  the  phase. 

In  the  implementation  of  this  block,  the  phase  ranges  from  -pi/2  to  3*pi/2,  instead  of  the  usual 
-pi  to  pi.  This  makes  it  possible  to  avoid  using  an  extra  IF  statement  in  the  program.  The  two 
forms  are  mathematically  equivalent. 

Inputs: 

i  In-Phase  component  of  input  data 

q  Quadrature  component  of  input  data 

Parameters: 

samples  Number  of  samples  in  input  and  output  vectors 


A-38 


Outputs: 

mag  Magnitude  of  input  data 

phase  Phase  of  input  data 


Name: 

histo/s_enonlinf 

Description: 

This  hierarchical  block  buffers  the  received  data  and  converts  it  from  serial  to  vector  data  and 
applies  it  to  the  histo/enonlin  nonlinearity.  In  addition,  the  summation  in  the  LO  decision 
statistic  is  formed  using  the  histo/correlator2  block.  The  output  of  the  correlator  is  converted 
back  to  serial  form. 

The  romelib/timing  block  generates  timing  waveforms  which  enable  and  disable  the 
nonlinearity  and  all  subsequent  blocks  at  appropriate  times.  Of  particular  interest  are  the 
CLK_OUT  and  HOLD_VEC  signals.  The  HOLD  VEC  signal  holds  vector  blocks  while  the 
serial  data  is  buffered  into  vector  form.  This  signal  goes  low  during  the  simulation  iteration  in 
which  the  vector  processing  completes  and  the  outputs  are  available.  The  CLK_OUT  signal 
provides  flexibility  for  use  in  multirate  systems  by  pulsing  high  when  a  valid  data  sample  is 
produced  at  the  output  of  the  histo/s_enonlin  block.  The  CLK_OUT  signal  has  the  same  period 
as  the  CLK_IN  signal,  thus  output  samples  are  produced  at  the  same  rate  as  that  of  the  input 
signal.  The  remaining  romelib/timing  signals  are  internal  to  the  histo/s_enonlin  block,  and  are 
used  in  serial  to  vector  and  vector  to  serial  buffering. 

The  SAMPLES_PER_SYMBOL,  NUM_REFS,  and  M_TYPE  are  characteristics  of  the 
reference  symbols.  The  reference  signals  must  be  created  before  using  this  block,  and  are  stored 
in  the  /spwdata/ref  library.  This  permits  the  designer  to  use  any  desired  signal  constellation. 
The  NUMJREFS  parameter  must  be  exported  to  the  histo/correlator2  block,  and  the  other  two 
parameters  are  exported  to  the  filename  of  the  rl/inf_vsource  block  as  "ref/"&substr 
("  rn_type":model,2,length("m_type":model)-2)&H_"&xstring(,,samples_per_symbor:model)& 
"spb_i’B  for  the  I  channel  signals  and  "ref/"&substr("m_type":model,2,length("m_type": 
model)-2)  &"_"&xstring(,'samples_per_symboi”:model)&,,spb_q’M  for  the  Q  channel  signals. 
For  example,  the  filenames  for  a  Binary  Phase  Shift  Keying  (BPSK)  signal  with  25  samples  per 
symbol  would  be  /spwdata/ ref/bpsk_25 spb_i  and  /spwdata/ref/bpsk_25spb_q. 

Inputs: 

in  Input  message  sequence  (complex  input) 

clkjn  Clock  input  for  romelib/timing 
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Parameters: 


samples 

Number  of  samples  in  input  vectors 

bins 

Number  of  bins  in  the  Probability  Density  Function  (PDF) 

samples  j5€r_symbol 

Number  of  samples  per  data  symbol 

nunwefs 

Number  of  reference  signals 

m_type 

Modulation  type 

Outputs: 

mnt 

Plot  of  the  Memoryless  Nonlinear  Transform 

g 

Complex  transformed  sequence 

out 

Decision  of  nonlinearity 

clk_out 

Timing  signal:  goes  high  when  valid  data  sample  is  available  at  the  output 
of  histo/s_enonlin 

hold_vec 

Timing  signal:  goes  low  when  vector  outputs  are  available 

pdf 

Plot  of  the  PDF  of  magnitude 

bp 

Breakpoints  (PDF  and  MNT  interval  boundaries) 

See  also: 

histo/enonlin,  histo/correlator2,  romelib/timing,  rl/inf_vsource 


Name: 

histo/s_enonlin_nct 

Description: 

This  hierarchical  block  buffers  the  received  data  and  converts  it  from  serial  to  vector  data  and 
applies  it  to  the  histo/enonlin  nonlinearity.  This  block  does  not  generate  the  summation  part  of 
the  LO  decision  statistic. 

The  romelib/timing  block  generates  timing  waveforms  which  enable  and  disable  the 
nonlinearity  and  all  subsequent  blocks  at  appropriate  times.  Of  particular  interest  are  the 
CLK_OUT  and  HOLD_VEC  signals.  The  HOLD_VEC  signal  holds  vector  blocks  while  the 
serial  data  is  buffered  into  vector  form.  This  signal  goes  low  during  the  simulation  iteration  in 
which  the  vector  processing  completes  and  the  outputs  are  available.  The  CLK_OUT  signal 
provides  flexibility  for  use  in  multirate  systems  by  pulsing  high  when  a  valid  data  sample  is 
produced  at  the  output  of  the  histo/s_enonlin_nc  block.  The  CLK_OUT  signal  has  the  same 
period  as  the  CLK_IN  signal,  thus  output  samples  are  produced  at  the  same  rate  as  that  of  the 
input  signal.  The  remaining  romelib/timing  signals  are  internal  to  the  histo/s_enonlin_nc  block, 
and  are  used  in  serial  to  vector  and  vector  to  serial  buffering. 
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The  SAMPLES_PER_SYMBOL,  NUM_REFS,  and  M_TYPE  are  characteristics  of  the 
reference  symbols.  The  reference  signals  must  be  created  before  using  this  block,  and  are  stored 
in  the  /spwdata/ref  library.  This  permits  the  designer  to  use  any  desired  signal  constellation. 
The  NUM_REFS  parameter  must  be  exported  to  the  histo/correlator2  block,  and  the  other  two 
parameters  are  exported  to  the  Filename  of  the  rl/inf_vsource  block  as  "ref/"&substr 
(',m_type":model,2,length("m_type":model)-2)&,,_"&xstring(,,samples_per_symbol":model)& 
"spb_i,M  for  the  I  channel  signals  and  " ref/ " &substr( " m_typeM : model ,2 , length( " m_typeM : 
model)-2)  &xstring( "  samples_per_sy mbol" :  model)& ”  spb_q’  “  for  the  Q  channel  signals. 

For  example,  die  filenames  for  a  Binary  Phase  Shift  Keying  (BPSK)  signal  with  25  samples  per 
symbol  would  be  /spwdata/ref/bpsk_25spb_i  and  /spwdata/ref/bpsk_25spb_q. 

Inputs: 

in  Input  message  sequence  (complex  input) 

clk_in  Clock  input  for  romelib/timing 

Parameters: 
samples 
bins 

sampies_per_symbol 

num_refs 
m_type 

Outputs: 

mnt  Plot  of  the  Memoryless  Nonlinear  Transform 

g  Complex  transformed  sequence  (vector) 

out  Complex  transformed  sequence  (serial) 

clk_out  Timing  signal:  goes  high  when  valid  data  sample  is  available  at  the  output 

of  histo/s_enonlin_nc 

hold_vec  Timing  signal:  goes  low  when  vector  outputs  are  available 

pdf  Plot  of  the  PDF  of  magnitude 

bp  Breakpoints  (PDF  and  MNT  interval  boundaries) 


See  also: 

histo/enonlin,  romelib/timing,  rl/inf_vsource,  histo/s_enonlin 


Number  of  samples  in  input  vectors 

Number  of  bins  in  the  Probability  Density  Function  (PDF) 

Number  of  samples  per  data  symbol 

Number  of  reference  signals 

Modulation  type 
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Name: 

histo/s_nonlin_ncf 

Description: 

This  hierarchical  block  buffers  the  received  data  and  converts  it  from  serial  to  vector  data  and 
applies  it  to  the  histo/nonlin6  nonlinearity.  This  block  does  not  generate  the  summation  part  of 
the  LO  decision  statistic. 


The  romelib/timing  block  generates  timing  waveforms  which  enable  and  disable  the 
nonlinearity  and  all  subsequent  blocks  at  appropriate  times.  Of  particular  interest  are  the 
CLK_OUT  and  HOLD_VEC  signals.  The  HOLD_VEC  signal  holds  vector  blocks  while  the 
serial  data  is  buffered  into  vector  form.  This  signal  goes  low  during  the  simulation  iteration  in 
which  the  vector  processing  completes  and  the  outputs  are  available.  The  CLK_OUT  signal 
provides  flexibility  for  use  in  multirate  systems  by  pulsing  high  when  a  valid  data  sample  is 
produced  at  the  output  of  the  histo/s_nonlin_nc  block.  The  CLK_OUT  signal  has  the  same 
period  as  the  CLK_IN  signal,  thus  output  samples  are  produced  at  the  same  rate  as  that  of  the 
input  signal.  The  remaining  romelib/timing  signals  are  internal  to  the  histo/s_nonlin_nc  block, 
and  are  used  in  serial  to  vector  and  vector  to  serial  buffering. 

The  SAMPLES_PER_SYMBOL,  NUM_REFS,  and  M_TYPE  are  characteristics  of  the 
reference  symbols.  The  reference  signals  must  be  created  before  using  this  block,  and  are  stored 
in  the  /spwdata/ref  library.  This  permits  the  designer  to  use  any  desired  signal  constellation. 
The  NUM_REFS  parameter  must  be  exported  to  the  histo/correlator2  block,  and  the  other  two 
parameters  are  exported  to  the  filename  of  the  rl/inf_vsource  block  as  "ref/"&substr 
( ”  m__type" :  model  ,2,length("  m_type" :  model)-2)&"_"  &xstring("  samples_per_symbol" :  model)& 
"spb_i’"  for  the  I  channel  signals  and  "ref/"&substr("m_type":model,2,length("m_type“: 
model)-2)  &"_"&xstring("samples_per_symbor:model)&"spb_q’"  for  the  Q  channel  signals. 
For  example,  the  filenames  for  a  Binary  Phase  Shift  Keying  (BPSK)  signal  with  25  samples  per 
symbol  would  be  / spwdata/ ref/bpsk_25 spb_i  and  /spwdata/ref/bpsk_25spb_q. 


Inputs: 

in  Input  message  sequence  (complex  input) 

clk_in  Clock  input  for  romelib/timing 


Parameters: 

samples 

bins 

samples_per_symbol 

num_refs 

mjype 


Number  of  samples  in  input  vectors 

Number  of  bins  in  the  Probability  Density  Function  (PDF) 

Number  of  samples  per  data  symbol 

Number  of  reference  signals 

Modulation  type 
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Outputs: 

mnt 

g 

out 

clk_out 

hold_vec 

pdf 

bp 


Plot  of  the  Memoryless  Nonlinear  Transform 
Complex  transformed  sequence  (vector) 

Complex  transformed  sequence  (serial) 

Timing  signal:  goes  high  when  valid  data  sample  is  available  at  the  output 
of  histo/s_nonlin_nc 

Timing  signal:  goes  low  when  vector  outputs  are  available 

Plot  of  the  PDF  of  magnitude 

Breakpoints  (PDF  and  MNT  interval  boundaries) 


See  also: 

histo/nonlinb,  romelib/timing,  rl/inf_vsource,  histo/serial_nonlin 


Name: 

histo/serial_nonlinf 

Description: 

This  hierarchical  block  buffers  the  received  data  and  converts  it  from  serial  to  vector  data  and 
applies  it  to  the  histo/nonlin6  nonlinearity.  In  addition,  the  summation  in  the  LO  decision 
statistic  is  formed  using  the  histo/correlator2  block.  The  output  of  the  correlator  is  converted 
back  to  serial  form. 

The  romelib/timing  block  generates  timing  waveforms  which  enable  and  disable  the 
nonlinearity  and  all  subsequent  blocks  at  appropriate  times.  Of  particular  interest  are  the 
CLK_OUT  and  HOLD_VEC  signals.  The  HOLD_VEC  signal  holds  vector  blocks  while  the 
serial  data  is  buffered  into  vector  form.  This  signal  goes  low  during  the  simulation  iteration  in 
which  the  vector  processing  completes  and  the  outputs  are  available.  The  CLK_OUT  signal 
provides  flexibility  for  use  in  multirate  systems  by  pulsing  high  when  a  valid  data  sample  is 
produced  at  the  output  of  the  histo/seriai_nonlin  block.  The  CLK_OUT  signal  has  the  same 
period  as  the  CLK_IN  signal,  thus  output  samples  are  produced  at  the  same  rate  as  that  of  the 
input  signal.  The  remaining  romelib/timing  signals  are  internal  to  the  histo/serial_nonlin  block, 
and  are  used  in  serial  to  vector  and  vector  to  serial  buffering. 

The  SAMPLES_PER_SYMBOL,  NUM_REFS,  and  M_TYPE  are  characteristics  of  the 
reference  symbols.  The  reference  signals  must  be  created  before  using  this  block,  and  are  stored 
in  the  /spwdata/ref  library.  This  permits  the  designer  to  use  any  desired  signal  constellation. 
The  NUM_REFS  parameter  must  be  exported  to  the  histo/correlator2  block,  and  the  other  two 
parameters  are  exported  to  the  filename  of  the  rl/inf_vsource  block  as  "ref/"&substr 
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( "  m_type  “ :  model  ,2 ,  length( "  m_type“ :  model)-2)& H  &xstring( "  samples_per_sy  mbol“ :  model)& 
"spbj’*  for  the  I  channel  signals  and  “  ref/ " &substr( “ m_type “ : model , 2 , length( “ m_type " : 
model)-2)  &"_"&xstring(MsamplesJ)er_symbo^:model)&,,spb_q,"  for  the  Q  channel  signals. 
For  example.  The  filenames  for  a  Binary  Phase  Shift  Keying  (BPSK)  signal  with  25  samples  per 
symbol  would  be  /spwdata/ref/bpsk_25spb_i  and  /spwdata/ref/bpsk_25spb_q. 

Inputs: 

in  Input  message  sequence  (complex  input) 

clk_in  Clock  input  for  romelib/timing 

Parameters: 
samples 
bins 

samples_per_  symbol 

num_refs 
mtype 

Outputs: 

mnt  Plot  of  the  Memoryless  Nonlinear  Transform 

g  Complex  transformed  sequence 

out  Decision  of  nonlinearity 

clk_out  Timing  signal:  goes  high  when  valid  data  sample  is  available  at  the  output 

of  histo/serial_nonlin 

hold_vec  Timing  signal:  goes  low  when  vector  outputs  are  available 

pdf  Plot  of  the  PDF  of  magnitude 

bp  Breakpoints  (PDF  and  MNT  interval  boundaries) 


See  also: 

histo/nonlin6,  histo/correlator2,  romelib/timing,  rI/inf_vsource 


Number  of  samples  in  input  vectors 

Number  of  bins  in  the  Probability  Density  Function  (PDF) 

Number  of  samples  per  data  symbol 

Number  of  reference  signals 

Modulation  type 
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Name: 

jam/channef 

Description: 

This  block  adds  Gaussian  noise  and  jammers  to  the  transmitted  data  sequence.  Currently  only 
three  Continuous  Wave  (CW)  jammers  and  a  Partial  Band  (PB)  are  implemented.  Any  other 
interference  blocks  may  be  easily  added.  (All  power  parameters  are  measured  in  dB) 

Inputs: 

in 

Transmitted  Message 

Parameters: 

J1_S 
phase  1 
freql 

J2_S 

phase2 

freq2 

J3_S 

phase3 

freq3 

Ratio  of  First  CW  jammer  power  to  signal  power 

First  CW  phase 

First  CW  frequency 

Ratio  of  Second  CW  jammer  power  to  signal  power 

Second  CW  phase 

Second  CW  frequency 

Ratio  of  Third  CW  jammer  power  to  signal  power 

Third  CW  phase 

Third  CW  frequency 

Jpb_S 

filt_order 

atten 

pfreq 

sfreq 

Ratio  of  PB  jammer  power  to  signal  power 

Lowpass  filter  order 

Passband  ripple  in  dB 

Passband  3  dB  edge  frequency 

Stopband  3  dB  edge  frequency 

Eb  No 

No 

Ratio  of  Bit  Energy  to  Gaussian  noise  power 

Gaussian  noise  power 

s_freq 

samples_per_symbol 

Sampling  frequency 

Number  of  Samples  per  Information  Symbol 

Outputs: 

out  Corrupted  Message 

See  also: 

rl/cw jammer,  rl/pbjammer 
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Name: 

linear/serial_lint 

Description: 

This  hierarchical  block  buffers  the  received  data  and  converts  it  from  serial  to  vector  data  and 
applies  it  to  the  histo/correlator2  block.  The  output  of  the  correlator  is  converted  back  to  serial 
form. 

The  romelib/ timing  block  generates  timing  waveforms  which  enable  and  disable  the 
nonlinearity  and  all  subsequent  blocks  at  appropriate  times.  Of  particular  interest  are  the 
CLK_OUT  and  HOLD_VEC  signals.  The  HOLD_VEC  signal  holds  vector  blocks  while  the 
serial  data  is  buffered  into  vector  form.  This  signal  goes  low  during  the  simulation  iteration  in 
which  the  vector  processing  completes  and  the  outputs  are  available.  The  CLK_OUT  signal 
provides  flexibility  for  use  in  multirate  systems  by  pulsing  high  when  a  valid  data  sample  is 
produced  at  the  output  of  the  linear/serial  Jin  block.  The  CLK_OUT  signal  has  the  same  period 
as  the  CLKJN  signal,  thus  output  samples  are  produced  at  the  same  rate  as  that  of  the  input 
signal.  The  remaining  romeiib/timing  signals  are  internal  to  the  linear/serial  Jin  block,  and  are 
used  in  serial  to  vector  and  vector  to  serial  buffering. 

The  SAMPLES_PER_SYMBOL,  NUM_REFS,  and  M_TYPE  are  characteristics  of  the 
reference  symbols.  The  reference  signals  must  be  created  before  using  this  block,  and  are  stored 
in  the  /spwdata/ref  library.  This  permits  the  designer  to  use  any  desired  signal  constellation. 
The  NUM_REFS  parameter  must  be  exported  to  the  histo/correlator2  block,  and  the  other  two 
parameters  are  exported  to  the  filename  of  the  rl/inf_vsource  block  as  "ref/‘*&substr 
( "  m  jype " :  model  ,2 ,  length("  m  Jype” :  model)-2)&B_M  &xstring( *  samples  j)er_symbo! " :  model)& 
"spbj’"  for  the  I  channel  signals  and  Href/"&substr("mjype":model,2,length("mjype'': 
modeI)-2)  &"_"&xstring("samples_per__symbor:model)«fe"spb_q’"  for  the  Q  channel  signals. 
For  example,  the  filenames  for  a  Binary  Phase  Shift  Keying  (BPSK)  signal  with  25  samples  per 
symbol  would  be  /spwdata/ref/bpsk_25spb_i  and  /spwdata/ref/bpsk_25spb_q. 

Inputs: 

in  Input  message  sequence  (complex  input) 

clkjn  Clock  input  for  romeiib/timing 

Parameters: 

samples 

samples_per_symbol 

num_refs 

mjype 


Number  of  samples  in  input  vectors 
Number  of  samples  per  data  symbol 
Number  of  reference  signals 
Modulation  type 
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Outputs: 

out 

clk_out 

hold_vec 

See  also: 
histo/comeiator2, 


Decision  of  linear  receiver 

Timing  signal:  goes  high  when  valid  data  sample  is  available  at  the  output 
of  linear/serialjin 

Timing  signal:  goes  low  when  vector  outputs  are  available 
romelib/timing,  rl/inf_vsource 
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Name: 

mipa/nonlin* 

Description: 

This  is  a  hierarchical  block  which  applies  the  nonlinearity  to  the  input  signal  in  vector  format. 
The  histo/rec_to_polar  block  converts  the  two  dimensional  input  signal  from  rectangular  to  polar 
coordinates,  and  a  MIPA  of  the  Probability  Density  Function  (PDF)  of  the  magnitude  is 
computed  by  the  mipa/pdf  block.  The  poly/mnt  block  obtains  the  Memoryless  Nonlinear 
Transform  (MNT)  from  the  MIPA  PDF  and  applies  it  to  the  magnitude  of  the  input  message 
sequence.  The  result  is  converted  back  to  rectangular  coordinates  by  the  histo/polar_to_rec 
block  using  the  unmodified  phase.  The  poly/plot  block  generates  a  plot  of  the  PDF  from  the 
polynomial  coefficients. 

Refer  to  the  poly/mnt  block  for  more  information  on  the  MNT. 


Inputs: 

ijn 

q_in 

In-Phase  component  of  input  message  sequence 
Quadrature  component  of  input  message  sequence 

Parameters: 

samples 

bins 

points 

Number  of  samples  in  input  vectors 

Number  of  bins  in  the  histogram  PDF 

Number  of  plot  points  for  PDF  and  MNT 

Outputs: 

pdf 

i_out 

q_out 

mnt 

bp 

CPA  PDF 

In-Phase  component  of  transformed  sequence 
Quadrature  component  of  transformed  sequence 
Plot  of  the  Memoryless  Nonlinear  Transform 
Breakpoints  (PDF  and  MNT  interval  boundaries) 

See  also: 

histo/rec_to_polar,  mipa/pdf,  poly/mnt,  poly/plot,  histo/polar_to_rec 
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Name: 

mipa/pdf 

Description: 

This  block  obtains  the  M-Interval  Polynomial  Approximation  (MIPA)  of  a  Probability  Density 
Function  (PDF).  The  MIPA  is  a  concatenation  of  polynomial  curves,  which  minimizes  the 
squared  error  between  the  approximation  and  the  actual  PDF.  This  implementation  only 
supports  Oth,  2nd,  and  4th  order  MIPAs. 


Inputs: 

message  Input  message  sequence 


Parameters: 

bins 

samples 

order 


Number  of  bins 

Number  of  samples  in  input  message  sequence 
Order  of  MIPA 


Outputs: 

bp  Breakpoints  (PDF  interval  boundaries) 

coeff  Polynomial  coefficients  values  ao,  a,,  ..,  aP 


Name: 

mipa/s_nonlin_nc+ 

Description: 

This  hierarchical  block  buffers  the  received  data  and  converts  it  from  serial  to  vector  data  and 
applies  it  to  the  mipa/nonlin  nonlinearity.  This  block  does  not  generate  the  summation  part  of 
the  LO  decision  statistic. 

The  romelib/timing  block  generates  timing  waveforms  which  enable  and  disable  the 
nonlinearity  and  all  subsequent  blocks  at  appropriate  times.  Of  particular  interest  are  the 
CLK_OUT  and  HOLD_VEC  signals.  The  HOLD_VEC  signal  holds  vector  blocks  while  the 
serial  data  is  buffered  into  vector  form.  This  signal  goes  low  during  the  simulation  iteration  in 
which  the  vector  processing  completes  and  the  outputs  are  available.  The  CLK_OUT  signal 
provides  flexibility  for  use  in  multirate  systems  by  pulsing  high  when  a  valid  data  sample  is 
produced  at  the  output  of  the  mipa/s_nonlin_nc  block.  The  CLK_OUT  signal  has  the  same 
period  as  the  CLK_IN  signal,  thus  output  samples  are  produced  at  the  same  rate  as  that  of  the 
input  signal.  The  remaining  romelib/timing  signals  are  internal  to  the  mipa/s_nonlin_nc  block, 
and  are  used  in  serial  to  vector  and  vector  to  serial  buffering. 
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The  SAMPLES_PER_SYMBOL,  NUM_REFS,  and  M_TYPE  are  characteristics  of  the 
reference  symbols.  The  reference  signals  must  be  created  before  using  this  block,  and  are  stored 
in  the  /spwdata/ref  library.  This  permits  the  designer  to  use  any  desired  signal  constellation. 
The  NUM_REFS  parameter  must  be  exported  to  the  histo/corre!ator2  block,  and  the  other  two 
parameters  are  exported  to  the  filename  of  the  rl/inf_vsource  block  as  "ref/“&substr 
("m_type“ :  mode!  ,2,length("  m_type''  :model)-2)&"_"&xstringCsamples_per_symbol“ :  model)& 
** spb_i ’ ’*  for  the  I  channel  signals  and  ,’ref/"&substr(“m_type“:model,2,length(“m_type“: 
model)-2)  &"_M&xstring("samples_per_symbor':model)&"spb_q’*  for  the  Q  channel  signals. 
For  example,  the  filenames  for  a  Binary  Phase  Shift  Keying  (BPSK)  signal  with  25  samples  per 
symbol  would  be  /spwdata/ref/bpsk_25spb_i  and  /spwdata/ref/bpsk_25spb_q. 

Inputs: 

in  Input  message  sequence  (complex  input) 

clkjn  Clock  input  for  romelib/timing 

Parameters: 
samples 
bins 

samples_per_symbol 

num_refs 
m_type 

Outputs: 

mnt  Plot  of  the  Memoryless  Nonlinear  Transform 

g  Complex  transformed  sequence  (vector) 

out  Complex  transformed  sequence  (serial) 

clk_out  Timing  signal:  goes  high  when  valid  data  sample  is  available  at  the  output 

of  mipa/s_nonlin_nc 

hold_vec  Timing  signal:  goes  low  when  vector  outputs  are  available 

pdf  Plot  of  the  PDF  of  magnitude 

bp  Breakpoints  (PDF  and  MNT  interval  boundaries) 


See  also: 

mipa/nonlin,  romelib/timing,  rl/inf_vsource,  mipa/serial_nonlin 


Number  of  samples  in  input  vectors 

Number  of  bins  in  the  Probability  Density  Function  (PDF) 

Number  of  samples  per  data  symbol 

Number  of  reference  signals 

Modulation  type 
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Name: 

mipa/serial_nonlint 

Description: 

This  hierarchical  block  buffers  the  received  data  and  converts  it  from  serial  to  vector  data  and 
applies  it  to  the  mipa/nonlin  nonlinearity.  In  addition,  the  summation  in  the  LO  decision  statistic 
is  formed  using  the  histo/correiator2  block.  The  output  of  the  correlator  is  converted  back  to 
serial  form. 


The  romelib/timing  block  generates  timing  waveforms  which  enable  and  disable  the 
nonlinearity  and  all  subsequent  blocks  at  appropriate  times.  Of  particular  interest  are  the 
CLK_OUT  and  HOLD_VEC  signals.  The  HOLD_VEC  signal  holds  vector  blocks  while  the 
serial  data  is  buffered  into  vector  form.  This  signal  goes  low  during  the  simulation  iteration  in 
which  the  vector  processing  completes  and  the  outputs  are  available.  The  CLK_OUT  signal 
provides  flexibility  for  use  in  multirate  systems  by  pulsing  high  when  a  valid  data  sample  is 
produced  at  the  output  of  the  mipa/serial_nonlin  block.  The  CLK_OUT  signal  has  the  same 
period  as  the  CLK_IN  signal,  thus  output  samples  are  produced  at  the  same  rate  as  that  of  the 
input  signal.  The  remaining  romelib/timing  signals  are  internal  to  the  mipa/serial_nonlin  block, 
and  are  used  in  seual  to  vector  and  vector  to  serial  buffering. 


The  SAMPLES  PER  SYMBOL,  NUM_REFS,  and  M_TYPE  are  characteristics  of  the 
reference  symbols.  The  reference  signals  must  be  created  before  using  this  block,  and  are  stored 
in  the  /spwdata/ref  library.  This  permits  the  designer  to  use  any  desired  signal  constellation. 
The  NUM_REFS  parameter  must  be  exported  to  the  histo/correlator2  block,  and  the  other  two 
parameters  are  exported  to  the  filename  of  the  rl/inf_vsource  block  as  Href/"&substr 
( "  m_type " :  model  ,2 ,  length( "  m_type " :  model)-2)& &xstring( "  samples_per_sy  mbol " :  model)& 
"spb_i’"  for  the  I  channel  signals  and  " ref/ " &substr( " m_type " : model , 2 , lengthf" m_type" : 
model)-2)  &"_"&xstring("samples_per_symbor:model)&"spb_q’"  for  the  Q  channel  signals. 
For  example,  the  filenames  for  a  Binary  Phase  Shift  Keying  (BPSK)  signal  with  25  samples  per 
symbol  would  be  /spwdata/ref/bpsk_25spb_i  and  /spwdata/ref/bpsk_25spb_q. 

Inputs: 

in  Input  message  sequence  (complex  input) 

clk_in  Clock  input  for  romelib/timing 


Parameters: 

samples 

bins 

samples_per_symbol 

num_refs 

m_type 


Number  of  samples  in  input  vectors 

Number  of  bins  in  the  Probability  Density  Function  (PDF) 

Number  of  samples  per  data  symbol 

Number  of  reference  signals 

Modulation  type 
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Outputs: 

mnt  Plot  of  the  Memoryless  Nonlinear  Transform 

g  Complex  transformed  sequence 

out  Decision  of  nonlinearity 

clk_out  Timing  signal:  goes  high  when  valid  data  sample  is  available  at  the  output 

of  histo/serial_nonlin 

hold_vec  Timing  signal:  goes  low  when  vector  outputs  are  available 

pdf  Plot  of  the  PDF  of  magnitude 

bp  Breakpoints  (PDF  and  MNT  interval  boundaries) 


See  also: 

mipa/nonlin,  histo/correlator2,  romelib/timing,  rl/inf_vsource 
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Name: 

poly/mnt* 

Description: 

This  block  obtains  the  Memory  less  Nonlinear  Transform  (MNT)  from  a  polynomial 
approximation  to  a  Probability  Density  Function  (PDF).  It  is  a  hierarchical  block  which  contains 
two  poly/mnt  calc  blocks.  One  of  them  computes  the  MNT  for  each  value  of  the  magnitude  R 
of  the  input  message  sequence,  and  the  other  generates  a  plot  of  the  MNT  from  the  minimum 
to  the  maximum  R  value  with  the  aid  of  a  vec/minmax_ramp  block.  This  block  only  operates 
on  the  magnitude  of  the  input  message  sequence.  For  a  P*  order  polynomial  the  MNT  is 

cfR1_  1  _  a,>2a2R>->PapRp-1 
R  a0+a1R+a2R2+--+apR p 

If  the  LT  parameter  is  set  to  1,  the  linear  transform  will  be  applied  to  the  coefficients.  That 
means  that  the  R  values  in  the  K*  bin  must  range  from  0  to  the  width  of  the  K*  bin 
(BP[K]-BP[K-1]),  whereas  without  the  transform  the  R  values  in  each  bin  will  range  from 
BP[K-1]  to  BP[K]. 


Inputs: 

bp 

coeff 

data 

Breakpoints  (PDF  and  MNT  interval  boundaries) 
PDF  polynomial  coefficients  values  a<,,  a„  aP 

Input  message  sequence  (1  or  2  dimensional) 

Parameters: 

samples 

bins 

points 

order 

It 

Number  of  samples  used  to  compute  MNT 

Number  of  bins 

Number  of  plot  points  in  G_PLOT 

Polynomial  order  (number  of  coefficients  -  1) 

Linear  transform  flag 

Outputs: 

g_val 

g_plot 

MNT  of  input  data 

Plot  of  MNT  from  minimum  to  maximum  data  value 

See  also: 

poly/mnt_calc,  vec/minmax_ramp 
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Name: 

poly/mnt_calc 

Description: 

This  function  computes  the  Memoryiess  Nonlinear  Transform  (MNT)from  a  polynomial 
approximation  to  a  Probability  Density  Function  (PDF)  for  each  value  of  the  magnitude  R  of  the 
input  message  sequence.  For  a  P*  order  polynomial  the  MNT  is 

1  aI+2a,R+™+PapRp-' 

gM--g-  — r  ---p 

R  ag+a^+a^+'+apR* 

In  the  case  of  a  one  dimensional  data  signal,  the  1/R  term  is  not  included  in  g[R].  For  a  two 
dimensional  data  signal,  this  block  operates  only  on  the  magnitude  R. 

If  the  LT  parameter  is  set  to  1 ,  the  linear  transform  will  be  applied  to  the  coefficients.  That 
means  that  the  R  values  in  the  K*  bin  must  range  from  0  to  the  width  of  the  K*  bin 
(BP[K]-BP[K-1]),  whereas  without  the  transform  the  R  values  in  each  bin  will  range  from 
BP[K-1]  to  BP[K], 

Inputs: 

bp  Breakpoints  (PDF  and  MNT  interval  boundaries) 

coeff  PDF  polynomial  coefficients  values  a<>,  at,  ...  ap 

data  Input  message  sequence  (1  or  2  dimensional) 

Parameters: 
bins 
samples 
order 
dim 
It 

Outputs: 

mnt  MNT  for  polynomial 


Number  of  bins 

Number  of  samples  used  to  compute  MNT 
Polynomial  order  (number  of  coefficients  -  1) 
Dimensions  of  data  signal  (1  or  2) 

Linear  transform  flag 
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Name: 

poly/mnt_calc_gt 

Description: 

This  function  computes  the  Memoryless  Nonlinear  Transform  (MNT)from  a  polynomial 
approximation  to  a  Probability  Density  Function  (PDF)  for  each  value  of  the  magnitude  R  of  the 
input  message  sequence.  For  a  P®  order  polynomial  the  MNT  is 

cfP1-  1  y2a,R+-..+Pa^-» 

R  a0+a1R+a2R2+-+apR p 

In  the  case  of  a  one  dimensional  data  signal,  the  1/R  term  is  not  included  in  g[R].  For  a  two 
dimensional  data  signal,  this  block  operates  only  on  the  magnitude  R. 

The  first  and  last  bins  are  fitted  with  the  MNT  of  a  Gaussian  PDF  instead  of  a  polynomial 
MNT.  THIS  BLOCK  IS  TO  BE  USED  IN  CONJUNCTION  WITH  THE  CPA/PDFjGT 
BLOCK. 

If  the  LT  parameter  is  set  to  1,  the  linear  transform  will  be  applied  to  the  coefficients.  That 
means  that  the  R  values  in  the  K®  bin  must  range  from  0  to  the  width  of  the  K*  bin 
(BP[K]-BP[K-1]),  whereas  without  the  transform  the  R  values  in  each  bin  will  range  from 
BP[K-1]  to  BP[K]. 

Inputs: 

bp  Breakpoints  (PDF  and  MNT  interval  boundaries) 

coeff  PDF  polynomial  coefficients  values  ao,  a„  ..,  aP 

data  Input  message  sequence  (1  or  2  dimensional) 

Parameters: 
bins 
samples 
order 
dim 
It 

Outputs: 

mnt  MNT  for  polynomial 


Number  of  bins 

Number  of  samples  used  to  compute  MNT 
Polynomial  order  (number  of  coefficients  -  1) 
Dimensions  of  data  signal  (1  or  2) 

Linear  transform  flag 
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Name: 

poly/mnt_gt+ 

Description: 

This  block  obtains  the  Memoryless  Nonlinear  Transform  (MNT)  from  a  polynomial 
approximation  to  a  Probability  Density  Function  (PDF).  It  is  a  hierarchical  block  which  contains 
two  poly/mnt_caic_gt  blocks.  One  of  them  computes  the  MNT  for  each  value  of  the  magnitude 
R  of  the  input  message  sequence,  and  the  other  generates  a  plot  of  the  MNT  from  the  minimum 
to  the  maximum  R  value  with  the  aid  of  a  vec/minmax_ramp  block.  This  block  only  operates 
on  the  magnitude  of  the  input  message  sequence.  For  a  P*  order  polynomial  the  MNT  is 

UP).  1  -  V2aiR*"-ParR'" 

R  a0+alR+a2R2+-+apRp 

The  first  and  last  bins  are  fitted  with  the  MNT  of  a  Gaussian  PDF  instead  of  a  polynomial 
MNT.  THIS  BLOCK  IS  TO  BE  USED  IN  CONJUNCTION  WITH  THE  CPA/PDF  GT 
BLOCK. 

If  the  LT  parameter  is  set  to  1,  the  linear  transform  will  be  applied  to  the  coefficients.  That 
means  that  the  R  values  in  the  K®  bin  must  range  from  0  to  the  width  of  the  K*  bin 
(BP[K]-BP[K-l])t  whereas  without  the  transform  the  R  values  in  each  bin  will  range  from 
BP[K-1]  to  BP[K]. 


Inputs: 

bp 

coeff 

data 

Parameters: 

bins 

samples 

order 

It 


Breakpoints  (PDF  and  MNT  interval  boundaries) 
PDF  polynomial  coefficients  values  a<j,  a„  ..,  aP 
Input  message  sequence  (1  or  2  dimensional) 


Number  of  bins 

Number  of  samples  used  to  compute  MNT 
Polynomial  order  (number  of  coefficients  -  1) 
Linear  transform  flag 


Outputs: 


g_val  MNT  of  input  data 

g_plot  Plot  of  MNT  from  minimum  to  maximum  data  value 

See  also: 

poly/mnt_calc_gt,  vec/minmaxjramp 


Name: 

poly/plot 

Description: 

This  function  plots  a  piecewise  polynomial  curve  with  the  number  of  intervals  equal  to  BINS. 
The  polynomial  in  each  interval  is  of  the  form 


y[X]  =  a0+a,X+a2X2+-+apxp 

where  P  is  the  ORDER.  The  X  value  is  varied  from  its  minimum  to  maximum  value  to  generate 
the  plot.  The  number  of  points  in  this  plot  is  equal  to  POINTS. 

If  the  LT  parameter  is  set  to  1,  the  linear  transform  will  be  applied  to  the  coefficients.  That 
means  that  the  X  values  in  the  K*  bin  must  range  from  0  to  (BP[K]-BP[K-1]),  whereas  without 
the  transform  the  X  values  in  each  bin  will  range  from  BP[K-1]  to  BP[K]. 

Inputs: 
bp 

coeff 

Parameters: 
bins 
points 
order 
It 

Outputs: 
plot 


X-axis  Breakpoints  (interval  boundaries) 
Polynomial  coefficients  values  ao,  a,,  ..,  aP 

Number  of  bins 

Number  of  points  to  plot 

Polynomial  order  (number  of  coefficients  -  1) 

Linear  transform  flag 


Polynomial  curve 


Name: 

poly/plot_gt 

Description: 

This  function  plots  a  piecewise  polynomial  curve  with  the  number  of  intervals  equal  to  BINS. 
The  polynomial  in  each  interval  is  of  the  form 


y[X]  =a0+a1X+a2X2+-+apx p 
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where  P  is  the  ORDER.  The  X  value  is  varied  from  its  minimum  to  maximum  value  to  generate 
the  plot.  The  number  of  points  in  this  plot  is  equal  to  POINTS. 

If  the  LT  parameter  is  set  to  1 ,  the  linear  transform  will  be  applied  to  the  coefficients.  That 
means  that  the  X  values  in  the  K*  bin  must  range  from  0  to  (BP[K]-BP[K-1]),  whereas  without 
the  transform  the  X  values  in  each  bin  will  range  from  BP[K-1]  to  BP[K]. 

The  first  and  last  bins  are  fitted  with  a  Gaussian  PDF  instead  of  a  polynomial. 

Inputs: 

bp  X  -axis  Breakpoints  (interval  boundaries) 

coeff  Polynomial  coefficients  values  a*,,  a,,  aP 

Parameters: 
bins 
points 
order 
It 

Outputs: 

plot  Polynomial  curve 


Number  of  bins 

Number  of  points  to  plot 

Polynomial  order  (number  of  coefficients  -  1) 

Linear  transform  flag 
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Name: 

rl/complex_cwt 

Description: 

This  block  generates  a  single  Continuous  Wave  (CW)  complex  sinusoid  with  a  specified 
amplitude,  frequency,  and  phase. 

Inputs: 


none 

Parameters: 

amp 

CW  amplitude 

phase 

CW  phase 

freq 

CW  frequency 

s_freq 

Sampling  frequency 

Outputs: 

cw_out 

Jammer  output 

Name: 

rl/complex_data+ 

Description: 

This  block  generates  a  random  Quadrature  Phase  Shift  Keying  (QPSK)  data  sequence.  A 
fundamental  relationship  exists  between  the  symbol  rate  and  the  sampling  frequency,  namely 
S  FREQ  =  Rs  *  SAMPLES_PER_SYMBOL.  This  means  that  all  three  of  these  parameters 
may  not  vary  independently,  but  one  must  be  a  function  of  the  other  two.  For  this  reason,  the 
symbol  rate  parameter  is  exported  as  Rs  =  S_FREQ  /  SAMPLES_PER_SYMBOL. 

Inputs: 

none 

Parameters: 

amp  Bit  amplitude 

Rs  Symbol  rate 

prob_zero  Probability  of  a  bit  being  zero 

s_freq  Sampling  frequency 
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Outputs: 

data_out  QPSK  data  stream 


Name: 

rl/cwjammer* 

Description: 

This  block  generates  three  Continuous  Wave  (CW)  jammers  with  specified  amplitude, 
frequency,  and  phase.  Any  number  of  CW  jammers  up  to  three  may  be  obtained  by  setting  the 
amplitude  of  the  undesired  jammers  to  an  arbitrarily  small  value. 

Inputs: 


none 

Parameters: 

ampl 

First  CW  amplitude 

phase  1 

First  CW  phase 

freql 

First  CW  frequency 

amp2 

Second  CW  amplitude 

phase2 

Second  CW  phase 

freq2 

Second  CW  frequency 

amp3 

Third  CW  amplitude 

phase3 

Third  CW  phase 

freq3 

Third  CW  frequency 

s_freq 

Sampling  frequency 

Outputs: 

cw_out 

Jammer  output 

See  also: 

rl/complex_cw 
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Name: 

rl/holdt 

Description: 

This  block  generates  an  output  that  goes  high  every  SAMPLES  iterations.  There  is  no  pulse 
on  the  Oth  iteration. 

Inputs: 

none 

Parameters: 

samples  Number  of  iterations  between  pulses 

Outputs: 

hold  Output  clock 


Name: 

rl/inf_vsourcef 

Description: 

This  block  is  identical  to  the  spb/vsource  SPW  library  block  except  that  when  an  end  of  file 
(EOF)  occurs  the  last  output  is  retained  for  the  remainder  of  the  simulation.  This  block  detects 
the  EOF  from  the  spb/vsource  and  holds  the  spb/vsource  from  then  on. 

Inputs: 

none 

Parameters: 
same  as  spb/vsource 

Outputs: 

same  as  spb/vsource 

See  also: 
spb/vsource 
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Name: 

rl/pbjammer1. 

Description: 

This  block  generates  a  baseband  Partial  Band  (PB)  jammer  by  low  pass  filtering  Gaussian 
noise  with  an  elliptic  lowpass  filter.  The  Passband  and  Stopband  edge  frequencies  are  assumed 
to  be  normalized  by  the  symbol  rate. 

Inputs: 


none 

Parameters: 

Jpb_S 

Ratio  of  jammer  power  to  signal  power  in  dB 

filt_order 

Lowpass  filter  order 

atten 

Passband  ripple  in  dB 

pfreq 

Passband  3  dB  edge  frequency 

sfreq 

Stopband  3  dB  edge  frequency 

s_freq 

Sampling  frequency 

Outputs: 

pb_out 

Jammer  output 

Name: 

rl/psk_err_cntt 

Description: 

THIS  BLOCK  DOES  NOT  WORK  CORRECTLY  BECAUSE  COMM/REAL  ERR  CNT 
DOES  NOT  WORK  CORRECTLY. 

Inputs: 

X  Actual  transmitted  message 

Y  Decision  of  receiver 

Parameters: 
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Outputs: 

pc  Probability  of  error 

num_symbois 

results_clk 

See  also: 
rl/real_err_cnt 


Name: 

rl/rcal_err_cntt 

Description: 

THIS  BLOCK  DOES  NOT  WORK  CORRECTLY  BECAUSE  COMM/REAL  ERR  CNT 
DOES  NOT  WORK  CORRECTLY. 

Inputs: 

X  Actual  transmitted  message 

Y  Decision  of  receiver 

Parameters: 


Outputs: 

pe  Probability  of  error 

num_symbols 

results_clk 

See  also: 

rl/real_err_cnt 
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Name: 

romelib/timingf 

Description: 

This  block  generates  the  timing  signals  for  the  nonlinear  receiver  and  all  subsequent  blocks. 
In  particular,  it  allows  for  serial-to-vector  conversion  of  the  input  data  to  the  nonlinearity  blocks, 
and  vector-to-serial  conversion  at  the  output.  The  CLK_IN  input  is  provided  to  facilitate  the 
integration  of  the  nonlinear  blocks  into  multirate  systems.  This  input  should  be  tied  to  a  timing 
source  which  goes  high  on  every  simulation  iteration  caring  which  a  valid  input  data  sample  is 
present.  The  only  parameter  for  this  block,  SAMPLES,  should  be  set  to  the  length  of  the  data 
vector  which  will  be  processed  by  the  nonlinear  block. 

The  sequence  of  output  signal  states  is  as  follows: 

a)  The  HOLD_IN  output  goes  low  every  iteration  that  the  CLK_IN  input  goes  high  and 
indicates  when  a  valid  input  signal  sample  is  present. 

b)  The  LOADOUT  output  goes  high  during  the  simulation  iteration  when  the  last  of  the 
SAMPLES  samples  is  available  at  the  input.  This  is  a  signal  to  the  circular_buffer  block 
that  all  the  required  data  samples  are  present  and  that  the  data  vector  is  complete. 

c)  The  HOLD_VEC  output  goes  low  one  simulation  iteration  after  the  last  of  the  SAMPLES 
samples  is  available  at  the  input.  This  signal  should  be  tied  to  the  HOLD  input  of  the 
nonlinearity  blocks  as  well  as  all  vector  output  sinks/blocks  associated  with  the 
nonlinearity  block,  allowing  these  blocks  to  go  active  during  the  simulation  iteration 
when  all  the  required  data  samples  are  present  in  vector  form. 

d)  The  LOADIN  goes  high  one  simulation  iteration  after  the  last  of  the  SAMPLES  samples 
is  available  at  the  input.  This  signal  is  used  to  "latch"  the  output  signal  vector  of  the 
nonlinearity  block  for  vector-to-serial  conversion. 

e)  The  HOLD_OUT  signal  goes  low  one  simulation  iteration  after  each  iteration  when  the 
CLK_IN  input  goes  high.  This  signal  is  provided  for  clocking  the  serial  blocks  located 
after  the  nonlinearity  block. 

f)  The  CLK_OUT  output  first  goes  high  on  the  same  simulation  iteration  as  the 
HOLD_VEC  output  goes  low  and  the  LOADIN  output  goes  high,  and  then  it  periodically 
goes  high  after  this  with  a  period  equal  to  that  of  the  CLK_IN  input  signal.  The 
CLK_OUT  signal  is  used  to  indicate  when  a  valid  serial  data  sample  is  present  at  the 
output  of  the  nonlinearity  block. 

Inputs: 

clkjn  Synchronizing  clock  input 

Parameters: 

samples  Length  of  input  data  vector 
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Outputs: 

hold_in  Activation  signal  for  serial  blocks  preceding  the  nonlinearity  block 

loadout  Indicates  when  input  data  vector  has  been  filled 

hold_vec  Activation  signal  for  the  vector  processing/output  blocks 

loadin  "Latches"  the  output  data  vector  for  vector-to-serial  conversion 

hold_out  Activation  signal  for  serial  blocks  following  the  nonlinearity  block 

clk_out  Synchronizing  clock  output 
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Name: 

vec/heap_sort 

Description: 

This  block  heap  sorts  the  input  data  in  ascending  order.  This  sort  is  of  order  N  log  (N).  For 
a  description  of  the  heap  sort  algorithm  see  Numerical  Recipes  in  C  (1988). 

Inputs: 

in  Input  vector 

Parameters: 

points  Number  of  points  in  the  vector 

Outputs: 

out  Sorted  vector 


Name: 

vec/minmax 

Description: 

This  block  finds  the  minimum  and  maximum  data  points  in  the  given  data  sequence. 
Inputs: 

in  Input  data  vector 

Parameters: 

length  Number  of  points  in  input  vector 

Outputs: 

min  Minimum  value 

max  Maximum  value 
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Name: 

vec/stretch_nl 

Description: 

This  block  plots  a  histogram  by  copying  the  value  of  each  bin  to  the  output  a  number  of  times 
that  is  proportional  to  the  width  of  the  bin.  The  total  width  of  the  histogram  is  T  WIDTH  = 
BPfBINS]  -  BP[0],  and  the  width  of  the  K*  bin  is  WIDTH [K]  =  BP[K]  -  BP[K-l].~The  K*  bin 
value  is  copied  to  the  output  N  times,  where  N  =  POINTS  *  WIDTHfK]  /  T_WIDTH. 

The  number  of  breakpoints  is  one  more  than  the  number  of  input  points. 

Inputs: 

bp  Breakpoints  of  histogram 

in  Input  histogram 

Parameters: 

bins  Number  of  bins  in  the  input  sequence 

points  Number  of  points  in  the  output  sequence 

Outputs: 

out  Plot  of  histogram 
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Name: 

vec/minmax_rampt 

Description: 

This  block  generates  a  linear  ramp  of  length  RAMP  _LEN  from  the  minimum  to  the  maximum 
values  in  the  given  data  sequence.  If  RAMP_LEN  =  1  it  is  not  possible  to  generate  a  line 
through  both  min  and  max  unless  min  =  max.  In  the  event  that  RAMP_LEN  =  1,  RAMP  is 
set  equal  to  the  average  of  min  and  max. 

Inputs: 
in 

Parameters: 
datalen 
ramp_len 

Outputs: 
ramp 

See  also: 

vec/minmax,  vec/ramp 


Input  data  vector 


Input  data  vector  length 
Ramp  vector  length 


Output  ramp 


Name: 

vec/ramp 


Description: 

This  block  generates  a  linear  ramp  of  length  POINTS  from  MIN  to  MAX.  If  POINTS  =  1 
it  is  not  possible  to  generate  a  line  through  both  MIN  and  MAX  unless  MIN  =  MAX.  In  the 
event  that  POINTS  =  1,  the  ramp  is  equal  to  the  average  of  MIN  and  MAX. 

Inputs: 

min  Minimum  data  value 

max  Maximum  data  value 


Parameters: 

points  Ramp  vector  length 

Outputs: 

ramp  Output  ramp 
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Name: 

vec/var 

Description: 

This  block  computes  the  mean  and  variance  of  the  input  data  vector  as  follows: 

.  LEN-l  .  LEN-1 

MEAN  =  _L-  £  IN[i]  VAR  =  — i— -  £  (IN[i] -MEAN)2 

LEN  to  LEN-1  " 

Inputs: 

in  Input  data  vector 

Parameters: 

len  Input  vector  length 

Outputs: 

mean  Mean  of  the  input  data  sequence 

var  Variance  of  the  input  data  sequence 
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Appendix  B 
SPW  Block  Diagrams 


This  Appendix  is  an  alphabetical  listing  of  the  details  of  the  hierarchical  IIT  SPW  blocks 
that  were  used  to  generate  the  results  in  Volume  1  of  this  report. 

When  the  BDE  prints  a  block  diagram,  the  parameters  do  not  reflect  their  exported 
values.  Rather,  the  default  values  for  these  parameters  are  displayed.  However,  in  the  actual 
simulation  the  parameter  values  have  all  been  correctly  exported. 
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ROMS  LABORATORY 


Rom*  Laboratory  plans  and  executes  an  interdiscipl inary 
program  in  research,  development,  test,  and  technology 
transition  in  support  of  Air  Force  Command,  Control, 
Communications  and  Intelligence  (C3I)  activities  for  all 
Air  Force  platforms.  It  also  executes  selected 
acquisition  programs  in  several  areas  of  expertise. 
Technical  and  engineering  support  within  areas  of 
competence  is  provided  to  ESC  Program  Offices  (POs)  and 
other  ESC  elements  to  perform  effective  acquisition  of 
C3I  systems.  In  addition,  Rome  Laboratory's  technology 
supports  other  AFMC  Product  Divisions,  the  Air  Force  user 
community,  and  other  DOD  and  non-DOD  agencies.  Rome 
Laboratory  maintains  technical  competence  and  research 
programs  in  areas  including,  but  not  limited  to, 
communications,  command  and  control,  battle  management, 
intelligence  information  processing,  computational 
sciences  and  software  producibility,  wide  area 
surveillance/sensors,  signal  processing,  solid  state 
sciences,  photonics,  electromagnetic  technology, 
superconductivity,  and  electronic 
reliability/maintainability  and  testability. 


